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ABSTRACT 

X 

A pump  system  was  built  and  instrumented  on  which  tests 
were  performed  yielding  single-  and  two-phase  actual  pump 
characteristics  for  first-quadrant  operation  (forward  flow 
and  forward  rotation)  of  two  different  centrifugal  impellers 
mounted  in  a simple  two-dimensional  volute.  Theoretical 
single-  and  two-phase  pump  characteristics  were  determined 
from  the  impeller  geometry,  the  Euler  equation,  and  the  use 
of  the  two-phase  flow  function.  The  head-loss  ratio,  the  loss 
of  head  in  two-phase  flow  divided  by  the  loss  of  head  in 
single-phase  flow,  was  then  plotted  versus  void  fraction. 

The  results  were  compared  to  an  earlier  theory  propose^ by 
J.  Mikielewicz  and  D.  G.  Wilsorc^whlcTi  predicted  the~head- 
loss  ratio  to  be  primarily  a function  of  void  fraction.  -Mjj, 
results  indicate  that  the  head-loss  ratio  is  a function  of 
flow  coefficient  as  well  as  void  fraction.  Flow-visualization 
studies  were  also  conducted  and  revealed  that  flow  regime 
inside  the  blade  passages  can  be  different  from  either  the 
inlet  or  outlet  flow  regimes,  and  greatly  affects  pump  per- 
formance. Further  experimentation,  on  a more  efficient  system 
employing  higher-specific-speed  impellers,  is  recommended  to 
study  more  closely  the  effect  of  flow  coefficient,  flow 
regime  and  void  fraction  upon  two-phase  pump  performance.  , 


Thesis  supervisor:  David  Gordon  Wilson 
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CHAPTER  1 


INTRODUCTION 


BACKGROUND 

r 

Safety  analyses  of  loss-of-coolant  accidents  (LOCA) 

in  pressurized-water  nuclear  reactors  (PWR)  require  the  pre- 
diction of  the  performance  of  the  main-coolant  pumps  in 
several  operating  modes  in  both  single-  and  two-phase  flow. 

In  a PWR,  sub-cooled  water  at  nominal  conditions  of  (typically) 
2250  psig  and  577  F is  circulated  through  the  reactor  core 
by  several  main-coolant  centrifugal  pumps,  each  in  its  own 
leg.  A LOCA  occurs  when  there  is  a break  in  one  or  more  of 
the  coolant  legs,  causing  the  sub-cooled  water  to  flash 
into  steam  while  flowing  through  the  coolant  pumps.  In 
addition,  the  flow  can  either  maintain  or  reverse  its  direct- 
ion, depending  upon  the  location  and  the  size  of  the  break. 
Thus,  prediction  of  the  behavior  of  centrifugal-pump  opera- 
tion in  two-phase  flow  is  necessary  before  the  thermal- 
hydraulic  response  of  a PWR  system  during  a LOCA  can  be 
fully  predicted. 

There  has  been  a continuing  effort  to  predict  more 
precisely  the  performance  of  centrifugal  pumps  during  two- 
phase  flow.  In  May  of  1964  D.  J.  Olson  of  the  Aerojet 
Nuclear  Company  (ANC)  published  the  results  of  testing  done 
on  a half-scale  centrifugal  pump,  which  was  part  of  a test 
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loop  designed  to  simulate  a LOCA  in  a PWR.  G.  L.  Sozzi  and 
G.  W.  Burnette  published  similar  test  results  in  November 
of  1971,  obtained  from  the  General  Electric  one-sixth-scale 
test  facility.  More  recently,  J.  Mikielewicz  and  D.  G. 
Wilson,  working  at  M.I.T.,  have  proposed  a method  of  predict- 
ing two-phase  pump  performance  as  a function  of  head-loss 
ratio  and  void  fraction.  This  method  has  been  applied  be 
T.  C.  Chan,  studying  at  M.I.T.,  to  the  ANC  data  with  reason- 
ably good  results. 

PURPOSE 

The  purpose  of  this  report  is  to  refine  the  method 
proposed  by  Mikielewicz  and  Wilson  by  investigating  whether 
the  head-loss  ratio  is  a function  of  other  variables  besides 
void  fraction.  In  addition,  flow-visualization  studies 
were  conducted  to  gain  an  insight  into  the  physical  phenom- 
enon of  the  flow  inside  the  blade  passages  and  volute,  and 
its  effects  on  pump  performance.  The  tests  and  conclusions 
contained  in  this  report  are  confined  to  first-quadrant 
operation  (forward  flow  and  forward  rotation). 


I 


15 


CHAPTER  2 

DESCRIPTION  OF  THE  APPARATUS 

GENERAL  DESCRIPTION 

The  tests  conducted  were  performed  on  an  air-water 
system.  City  water  at  a pressure  of  50  ± 5 psig  was 
supplied  by  a 1.61-inch  I.D.  pipe  while  the  air  was  supplied 
by  the  lab-air  supply  line  through  a 0.25-inch  nylon  fitting 
in  the  inlet  piping.  The  test  facility  employed  two  differ- 
ent impellers,  mounted  vertically  in  a clear  Plexiglass 
casing,  and  powered  by  a variable-speed  A.C.  motor.  The 
casing  emptied  into  the  botoom  of  a 55-gallon  drum  which  was 
fitted  with  a metal  weir  to  allow  for  constant  backpressure 
during  operation,  and  having  a 4-inch-diameter  drain  on  the 
opposite  side  of  the  weir.  A venturi  was  used  to  measure 
water-flow  rate  while  the  air  flow  was  measured  by  rotameters. 
Static  head  across  the  impeller  was  measured  by  a water 
manometer  and  supply  air  pressure  was  measured  by  a standard 
pressure  gauge.  An  overall  diagram  of  the  apparatus  is 
given  in  figure  2-1  and  a picture  of  the  actual  test  rig  is 
given  in  photo  2-1. 

IMPELLERS 

The  first  impeller  used  in  the  test  rig  was  a low- 
specific-speed,  low-head,  molded-plastic  centrifugal  impeller 
with  backward-swept  blades.  It  had  twenty-four  blades  and  a 


1G 


WATFR  SUPPLY  LINE  

I/?"  GLOBE  VALVE  — 
rj‘,  GAL  I ON  DRUM  — 


- VENTURI 

1/2"  GLOBE  VALVE 


A"  DRAIN 


22"  METAL 
WEIR 


ROTAMETER 


FIGURE  2-1.  OVERALL  VIEW  OF  THE  M.I.T.  FLOW-VISUALIZATION 
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blade  outlet  angle  of  46°.  The  blades  were  0.5  inches 
thick  at  the  tips  and  maintained  a constant  height  of  0.5 
inches  from  shroud  to  tip.  The  eye  of  the  impeller  was 
4.25  inches  in  diameter  while  the  blade-tip  diameter  was 
8.19  inches  (see  photos  2-2  and  2-3).  The  plastic  impeller 
was  mounted  on  an  aluminum  backplate  to  prevent  warping  at 
high  speed. 

The  second  impeller  used  was  a low-specific-speed, 
shrouded,  brass  centrifugal  impeller  with  backward- swept 
blades.  This  impeller  had  five  blades  and  a blade  outlet 
angle  of  25°.  The  blades  were  0.25  inches  thick  and  0.385 
inches  high  at  the  tips.  The  impeller  eye  was  2.81  inches 
in  diameter  while  the  blade-tip  diameter  was  7.31  inches. 

A detailed  drawing  of  the  impeller  is  given  in  figure  2-2. 
CASINGS 

The  original  impeller  casing  consisted  of  three  30-inch 
square  pieces  of  Plexiglass:  a 0.5-inch  thick  scroll  piece 
encased  between  a 0.5-inch  thick  backplate  and  a 0.25-inch 
thick  coverplate.  The  Plexiglass  casing  allowed  visual  obser- 
vation of  the  flow  both  within  the  impeller  passages  and  in 
the  scroll. 

The  scroll  was  designed  for  optimum  first-quadrant 
operation  of  the  original  impeller.  This  was  done  by  first 
assuming  a best  efficiency  value  of  <f>  for  the  impeller 
(reference  2-1).  Based  upon  this  value  of  <J>,  the  impeller 
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geometry  and  a slip-angle  calculated  from  the  Busemann  slip- 
factor  correlation  (reference  2-2),  the  0.5-inch  deep  scroll 
increased  0.388  inches  in  radius  for  every  15°  of  arc  to  a 
total  outlet  area  of  4.66  square  inches.  This  scroll  config- 
uration allowed  the  absolute  velocity  of  the  flow  at  the 
blade  tips  to  be  maintained  throughout  the  scroll.  The 
details  and  calculations  involved  in  determining  the  scroll 
shape  are  contained  in  Appendix  B. 

The  inlet  piping  to  the  casing  consisted  of  a vertical 
1. 5-inch-I.D.  Plexiglass  tube  joined  to  the  coverplate.  The 
casing  and  scroll  emptied  into  the  base  of  a 55-gallon  drum. 

An  overall  view  of  the  scroll  configuration  is  given  in 
figure  2-3. 

The  casing  used  with  the  second  impeller  consisted  of 
the  original  backplate  and  scroll.  However,  two  0.75-inch 
thick  coverplates  were  used  to  seal  the  impeller  more  securely 
in  the  casing.  A detailed  cross-sectional  drawing  of  the 
second  impeller  in  the  casing  is  given  in  figure  2-4. 

The  inlet  piping  was  increased  to  2 . 5-inches-I . D . Plex- 
iglass tubing  to  match  more  closely  the  eye  of  the  new  im- 
peller, while  the  scroll  exit  conf iguration  remained  the  same. 
PIPING  CONFIGURATION 

Figure  2-1  and  photo  2-1  show  the  overal  system  piping 
configuration  for  the  plastic  impeller  tests.  The  water  was 
supplied  by  a 1 . 61-inch-I .D.  pipe  at  50  ± 5 psig.  The  water 


FIGURE  2-3 
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flow  was  controlled  by  means  of  two  1.5-inch  globe  valves. 

The  first  valve,  located  upstream  of  the  venturi,  was  reg- 
ulated to  set  the  majority  of  the  flow  rates,  while  the 
valve  downstream  of  the  venturi  remained  fully  open,  except 
when  very  low  water-flow  rates  of  approximately  10  gallons 
per  minute  were  desired.  In  such  cases,  both  valves  were 
used  to  regulate  the  water  flow  so  that  the  water  manometer 
connected  to  the  venturi  could  be  read  more  easily. 

The  venturi  was  located  in  a 16.5-inch-long  horizontal 
section  of  1 . 61-inch-I . D . pipe  with  the  entrance  plane  9.47 
inches  from  the  forward  elbow,  and  the  exit  plane  2.7  inches 
from  the  rear  elbow.  The  venturi  was  inserted  into  the  pipe 
and  held  in  place  by  epoxy  cement  and  the  valve  connection 
at  the  throat.  The  pressure  taps  were  located  3.25  inches 
ahead  of  the  entrance  plane  and  at  the  mid-point  of  the 
throat.  0 . 25-inch-O.D.  Polyflo  tubing  connected  the  1/4- 
inch  globe  vlaves  at  the  pressure  taps  to  the  60-inch  water 
manometer . 

A 12-inch-long  rubber  hose  connected  the  steel  supply 
piping  to  the  10-inch-high,  1 . 5-inch-I . D.  Plexiglass  inlet 
pipe.  This  rubber  hose  also  contained  a 5-inch  straw-bundle 
flow  straightener  which  was  used  to  minimize  inlet  swirl. 

The  Plexiglass  inlet  pipe  had  a 1/32-inch  pressure  tap 
located  3 inches  upstream  of  the  impeller  inlet  to  measure 
inlet  static  head.  This  pressure  tap  was  connected  by  means 
of  0. 25-inch-O.D.  Polyflo  tubing  to  a 30-inch  water  manometer. 
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The  pump-outlet  static-head  pressure  tap  was  located  at 
mid-stream  in  the  exit  volute,  7 inches  from  the  entrance 
to  the  55-gallon  drum,  and  similarly  connected  to  the  30- 
inch  manometer. 

The  air-supply  system  consisted  of  the  1/2-inch-I . D. 
air-supply  pipe  being  reduced  to  1/4-inch,  being  fed  through 
a rotameter  and  pressure  gauge,  and  finally  injecting  air 
into  the  inlet  steel  piping  through  a 1/8-inch  Polyflo  fit- 
ting. The  air-inlet  fitting  was  located  22  inches  upstream 
of  the  impeller  eye. 

The  piping  configuration  remained  basically  the  same 
for  the  bronze  impeller  tests.  The  only  alteration  was  that 
the  Plexiglass  inlet  pipe  now  had  an  inside  diameter  of 
2.5  inches. 

INSTRUMENTATION 

The  water-flow  rate  was  measured  by  a venturi  connected 
to  a 60-inch  water  manometer.  The  venturi  was  designed 
(reference  2-3)  and  built  to  accomodate  a flow  range  of  10 
to  70  gallons  per  minute.  The  throat  diameter  was  1.258 
inches  and  the  entrance  and  exit  planes  had  a diameter  of 
1.544  inches.  The  inlet  angle  measured  5.52°  while  the 
outlet  angle  was  3.98°.  The  venturi  was  made  from  machined 
aluminum  stock.  Complete  details  and  dimensions  of  the 
venturi  are  given  in  figure  2-5. 

The  air-flow  rate  was  measured  by  two  separate 
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rotameters,  The  smaller  rotameter,  used  for  flow  rates  of 
0.002  to  0.014  cubic  feet  per  second,  was  a Fischer-Porter 
1/4  - 20  - G - 5/81  glass  tube  with  a black-bead  float. 

The  larger  rotameter  covered  a flow  range  of  0.02  to  0.10 
cubic  feet  per  second,  and  was  a Fischer-Porter  1/2  - 27  - 
G - 10/80  glass  tube  with  a 0.5-inch  SVT-45  stainless-steel 
float.  Through  the  use  of  the  two  rotameters  and  the  venturi, 
void  fractions  from  0.0  to  0.62  were  obtained.  The  details 
of  the  void-fraction  calculation  are  discussed  in  Chapter  3. 

The  inlet  air  pressure  was  measured  by  an  Acco  Heli- 
coid air-pressure  gauge.  The  gauge  had  a range  of  0.0  to 
60.0  psig  and  was  used  to  record  inlet  air  pressure  necessary 
to  determine  the  slip-ratio,  s,  and  the  two-phase  density, 

Ptp* 

The  static  head  across  the  impeller  was  measured  by 
means  of  a 30-inch  water  manometer.  This  manometer  was 
connected  by  0. 25-inch-O.D.  Polyflo  tubing  to  1/32-inch 
pressure  taps  located  in  the  Plexiglass  inlet  pipe  and  the 
exit  scroll  channel.  Thus,  the  static-head  rise  across  the 
impeller  was  read  directly  from  the  30-inch  manometer  in 
inches  of  water. 

Finally,  the  angular  velocity  of  the  impeller  was 
determined  by  matching  its  rate  of  rotation,  controlled  by 
a variable-speed  electric  motor,  with  a setting  on  a Strobotac. 
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CALIBRATION 

The  calibration  of  the  venturi  was  conducted  while  it 
was  physically  in  the  test  loop,  although  a slight  modifica- 
tion of  the  loop  was  necessary.  The  short  rubber-hose  con- 
nection containing  the  flow  straightener  and  connecting  the 
metal  pipe  to  the  Plexiglass  inlet  pipe  was  removed.  A long 
rubber  hose  was  connected  to  the  piping  immediately  follow- 
ing the  second  globe  valve,  and  led  to  a sump  at  atmospheric 
pressure.  The  water  flow  was  then  turned  on  and  set  to  a 
steady-state  value  on  the  venturi  manometer  while  exhausting 
into  the  sump.  Once  steady  flow  was  achieved  the  large 
rubber  hose  was  swung  into  a 55-gallon  catch  tank  and  the  time 
to  fill  the  drum  was  recorded  by  a stop  watch.  Numerous 
values  of  AH  were  set  on  the  venturi  manometer  and  the  time 
to  fill  the  drum  was  recorded.  From  these  data  points  a 
calibration  curve  for  the  venturi  was  constructed.  This 
curve  is  shown  in  figure  2-6. 

The  two  separate  rotameters  used  during  the  testing 
were  also  calibrated  while  physically  in  the  test  loop. 

Each  rotameter  was  calibrated  separately  using  the  same 
method.  A five-gallon  plastic  container  was  filled  with  water 
and  submerged  upside  down  in  a water  sump.  The  air  flow 
was  turned  on  and  a steady  reading  was  set  on  the  rotameter 
in  the  test  loop.  Tubing,  0.25-inch  Polyflo,  was  run  from 
the  test  loop  behind  the  pressure  gauge  and  inserted  inside 


METER  AH,  (in.  of  H20) 


Ql,  FLOW  RATE, 


gal. /min. 
(ft.Vsec.) 


FIGURE  2-6 


CALIBRATION  CURVE  FOR  THE  VENTURI  AND 
WATER— FLOW- RATE  MANOMETER 
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the  submerged  container.  The  time  taken  to  displace  the 
known  volume  of  air  was  recorded  for  numerous  settings  on 
the  rotameters,  and  calibration  curves  were  constructed  for 
both  rotameters  (figure  2-7  and  figure  2-8) . 

The  Acco  Helicoid  air-pressure  gauge  was  calibrated 
on  an  Ashcroft  portable  dead-weight  tester.  The  gauge  read- 
ing was  found  to  be  constantly  0.5  psig  below  the  true  value 
throughout  the  range  of  the  gauge. 

AUXILIARY  EQUIPMENT 

The  pump  impeller  was  powered  by  a Master  variable- 
speed  A . C . electric  motor.  The  motor  had  a speed  range  of 
180  to  1725  revolutions  per  minute.  The  motor  shaft  was 
connected  to  the  impeller  shaft  at  a right  angle  by  means 
of  two  bevel  gears  with  tooth  angles  of  23°.  The  motor  was 
bolted  to  a wooden  baseplate  to  minimize  vibrations.  The 
impeller  shaft  was  held  at  the  impeller  end  by  a 1 3/16-inch 
stainless-steel  ball  bearing,  while  at  the  opposite  end  it 
was  connected  to  a 1/2-inch  bearing  attached  to  the  wooden 
baseplate. 

The  55-gallon  drum  into  which  the  scroll  emptied  was 
used  to  provide  a constant  back  pressure  during  operation, 
through  the  use  of  the  22-inch-high  metal  weir.  Flow  from 
the  scroll  exit  entered  the  drum  on  one  side  of  the  weir  and 
built  up  in  height  until  spillage  over  and  around  the  weir 
equalled  the  input  flow.  The  spillage  exited  the  drum  on  the 
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Qv,  AIR-FLOW  RATE  (ft3/sec) 

FIGURE  2-7.  CALIBRATION  CURVE  FOR  THE  FISCHER  & PORTER 
1/4  - 20  - G - 5/Cl  ROTAMETER  WITH  BLACK- 
GLASS-BEAD  FLOAT 


ROTAMETER  READING 


FIGURE  2-8.  CALIBRATION  CURVE  FOR  TIIE  FISCIIER  & PORTER 
1/2  - 27  - G - 10/C0  ROTAMETER  WITH 
1/2  - SVT  - 45  STAINLESS-STEEL  FLOAT 
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other  side  of  the  weir  through  the  4-inch  drain  hole.  For  a 
constant  input  flow  the  height  of  the  water  in  the  drum 
remained  constant,  thus  providing  a constant  back  pressure. 

A constant  back  pressure  was  necessary  to  take  accurate  read- 
ings both  across  the  impeller  and  for  the  venturi. 
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CHAPTER  3 

EXPERIMENTAL  PROCEDURE 

SINGLE-PHASE  CHARACTERISTIC  TESTS 

With  the  plastic  impeller  in  use,  26  single-phase 
data  points  were  obtained  while  28  points  were  obtained  us- 
ing the  bronze  impeller.  The  method  used  to  obtain  the  data 
from  both  of  the  impellers  was  the  s .«te.  A constant  speed 
was  set  on  the  variable-speed  motor  and  checked  by  the  use 
of  the  Strobotac.  The  water-flow  rate  was  then  adjusted  to 
an  appropriate  value  on  the  60-inch  manometer,  and  determined 
from  the  venturi  calibration  curve.  For  each  setting  of 
pump  speed  and  water- flow  rate  the  corresponding  static 
head  across  the  pump  was  read  off  the  30-inch  manometer. 
Speeds  of  190,  250,  318  and  390  rev/min  and  flow  rates  of 
0.0  to  46.0  gal/min  were  tested  using  the  plastic  impeller. 
With  the  bronze  impeller,  speeds  of  200,  300,  350  and  400 
rev/min  and  flow  rates  of  0.0  to  50.0  gal/min  were  tested. 
Care  was  taken  to  allow  the  back  pressure  to  reach  a steady 
value  (constant  water  height  in  the  exit  drum)  and  to  elimin- 
ate any  air  bubbles  in  the  manometer  lines  before  any  read- 
ings were  taken. 

The  shut-off  head  was  determined  for  the  various  speeds 
by  blocking  the  drum  drain  hole,  shutting  off  the  flow,  and 
allowing  the  water  in  the  drum  and  in  the  inlet  piping  to 
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achieve  a level  configuration.  The  motor  was  then  turned  on, 
pumping  the  water  to  a higher  level  in  the  drum  and  lowering 
the  level  in  the  inlet  piping.  When  a steady-state  level 
was  attained,  the  difference  in  height  between  the  water 
level  in  the  drum  and  that  viewed  in  the  Plexiglass  inlet 
pipe  was  the  shut-off  head. 

Since  the  correlation  proposed  by  Mikielewicz  and 
Wilson  is  based  upon  a head- loss  ratio  employing  total  heads, 
the  dynamic  head  across  the  pump  had  to  be  calculated  and 
added  to  the  measured  static  heads  in  order  to  obtain  the 
actual  total  single-phase  head,  AHQsp  . The  dynamic  head  was 
calculated  in  the  following  manner. 


H 


dyn  sp 


where  H 


dyn  sp 


2 2 

= cf  . (Q)  x L_ 

2g  A 2g 

= single-phase  dynamic  head. 


Therefore, 
H 


dyn  spl 


= <2-)  - i- 

'a  ' 5? 


H _ ,Q  X x 1 

"dyn  sP2  2g 


and 


AH 


dyn  sp 


Hdyn  sp2  Hdyn  spl  *5g^  ^ 2 “ 2 

A2  A1 


] 


where  for  the  plastic  impeller  Aj  = 0.0123  fti 


A , = 0.0324  ft 
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and  for  the  bronze  impeller  A^  = 0.0341  ft 

A2  = 0.0324  ft2 

Having  obtained  a value  of  AH  for  each  of  the  data 

osp 

points  recorded,  a plot  of  head  coefficient,  ip' , versus  flow 
coefficient,  <p , was  then  made. 

gAH0 

Head  coefficient  = i|>'  = — =— 

u 

c 

Flow  coefficient  = <b  = — 

u 

The  first-quadrant  single-phase  pump  characteristic  curves 
obtained  by  actual  data  plotting  are  shown  in  figures  3-1  and 
3-2  for  the  plastic  and  bronze  impellers  respectively.  The 
single-phase  data  used  to  construct  these  curves  are  listed 
in  tables  3-1  and  3-2  respectively. 

The  data  from  the  bronze  impeller  were  run  on  a curve- 
fitting computer  program  (reference  3-1,  Appendix  E)  and  the 
resulting  curve- fitted  plot  is  shown  in  figure  3-3.  The 
equation  for  the  curve  in  figure  3-3  is  as  follows. 

ip ' sp  = 0.413  + 0.580  <p2  - 5.772  4>22  - 614.27  <P2 3 + 

17931.1  <f>24  - 48526.0  4> 2 5 + 106810.0  <j>26  (3-1) 

TWO-PHASE  CHARACTERISTIC  TESTS 

Before  two-phase  data  could  be  obtained  from  the  test 
facility  a method  of  calculating  void  fraction  had  to  be 
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SINGLE-PHASE  CHARACTERISTIC 
OF  PLASTIC  PUMP 
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FIGURE  3-1. 
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FIGURE  3-2.  SINGLE-PHASE  CHARACTERISTIC  OF  THE  SRONSE 
IMPELLER 


39 


TABLE  3-1 


SINGLE-PHASE  PLASTIC- IMPELLER  DATA 


MOTOR 

SPEED 

(rpm) 

TEST 

NO. 

WATER- 
FLOW  RATE 

( f t^/sec) 

TOTAL 

HEAD 

(in.  II20) 

FLOW 

COEFF. 

4> 

HEAD 

COEFF. 

r 

190 

01 

0.0 

6.25 

0.0 

0.3640 

02 

0.0223 

3.78 

0.0787 

0.220 

03 

0.0361 

2.43 

0.1276 

0.1412 

04 

0.0512 

0.03 

0.1811 

0.0017 

250 

05 

0.0 

11.0 

0.0 

0.370 

06 

0.023 

7.58 

0.0598 

0.2549 

07 

0.0361 

5.63 

0.0969 

0.1892 

08 

0.0512 

3.23 

0.1376 

0.1086 

09 

0.0624 

1.39 

0.1676 

0.0468 

318 

10 

0.0 

17.75 

0.0 

0.3690 

11 

0.0223 

13.38 

0.0470 

0.2781 

12 

0.0361 

11.23 

0.0762 

0.2334 

13 

0.0512 

8.73 

0.1082 

0.1815 

14 

0.0624 

6.09 

0.1317 

0.1267 

15 

0.0719 

4.17 

0.1515 

0.0868 

16 

0.0806 

1.55 

0.1703 

0.0323 

390 

17 

0.0 

26.0 

0.0 

0.3590 

18 

0.0223 

21.73 

0.0384 

0.0300 

19 

0.0361 

18.33 

0.0621 

0.2533 

20 

0.0512 

15.93 

0.0882 

0.2202 
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MOTOR 

SPEED 

(rpm) 

390 


TABLE  3-1 

(continued) 

TEST 

NO. 

WATER- 
FLOW  RATE 

( ft^/sec) 

TOTAL 

HEAD 

(in.  H20) 

FLOW 

COEFF. 

HEAD 

COEFF. 

r 

21 

0.0624 

13.39 

0.1074 

0.1851 

22 

0.0719 

10.77 

0.1239 

0.1489 

23 

0.0806 

8.74 

0.1389 

0.1208 

24 

0.0888 

6.20 

0.1527 

0.0857 

25 

0.0950 

4.48 

0.1627 

0.0620 

26 

0.1024 

2.07 

0.1761 

0.0285 

1 
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TABLE  3-2 


SINGLE-PHASE 

BRONZE- IMPELLER 

DATA 

MOTOR 

SPEED 

(rpm) 

TEST 

NO. 

WATER- 
FLOW  RATE 

( f t Vsec) 

TOTAL 

HEAD 

(in.  H20) 

FLOW 

COEFF. 

<P 

HEAD 

COEFF. 

200 

01 

0.0 

6.5 

0.0 

0.4283 

02 

0.0361 

4.32 

0.0974 

0.2847 

03 

U . 0512 

0.25 

0.1381 

0.0165 

04 

0.0624 

-3.73 

0.1684 

-0.2458 

05 

0.0719 

-9.11 

0.1940 

-0.6003 

300 

06 

0.0 

14.13 

0.0 

0.4138 

07 

0.0361 

12.0 

0.0649 

0.3514 

08 

0.0512 

9.4 

0.0921 

0.2753 

09 

0.0624 

5.3 

0.1122 

0.1552 

10 

0.0719 

0.39 

0.1293 

0.0114 

11 

0.0806 

-4.39 

0.1450 

-0.1286 

12 

0.0886 

-9.57 

0.1594 

-0.2803 

350 

13 

0.0361 

16.52 

0.0557 

0.3555 

14 

0.0512 

14.05 

0.0789 

0.3023 

15 

0.0624 

11.37 

0.0962 

0.2446 

16 

0.0719 

7.19 

0.1109 

0.1547 

17 

0.0806 

3.11 

0.1243 

0.0669 

18 

0.0886 

-1.67 

0.1366 

-0.0359 

19 


0.0953 


-6.55 


0.1469 


-0.1409 
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TABLE  3-2 

(continued) 

ilOTOR 

SPEED 

(rpm) 

TEST 

NO. 

WATER- 
FLOW  RATE 

( ftVsec) 

TOTAL 

HEAD 

(in.  H20) 

FLOW 

COEFF. 

HEAD 

COEFF. 

400 

20 

0.0 

24.63 

0.0 

0.4057 

21 

0.0361 

22.62 

0.0487 

0.3726 

22 

0.0512 

20.35 

0.0691 

0.3352 

23 

0.0624 

17.17 

0.0842 

0.2829 

24 

0.0719 

13.29 

0.0970 

0.2189 

25 

0.0806 

9.51 

0.1087 

0.1567 

26 

0.0886 

5.03 

0.1195 

0.0829 

27 

0.0953 

0.45 

0.1286 

0.0074 

28 

0.1022 

-4.42 

0.1379 

-0.0728 

CURVE  — FITTED  HEAD  COEFFICIENT 


FLOW  COEFFICIENT 
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established.  Since  the  air  and  water  may  travel  at  differ- 
ent velocities  throughout  the  system,  the  void  fraction  can- 
not be  inferred  from  a ratio  of  the  respective  inlet  volume 
flows.  The  method  chosen  for  calculating  the  void  fraction 
during  these  tests  was  based  on  the  drift-flux  model  of  two- 
phase  flow  developed  by  Zuber  and  Wallis  (reference  3-2) . 

The  drift-flux  model  satisfactorily  accounts  for  the  influ- 
ence of  mass  velocity  on  the  void  fraction  as  seen  in  the 
separated-f low  model,  and  is  useful  in  the  bubbly-,  slug-, 
and  churn-flow  regimes.  These  flow  regimes  include  void 
fractions  from  0.0  to  approximately  0.80.  The  specific 
equation  developed  from  the  drift-flux  model  used  to  calcu- 
late void  fractions  throughout  the  experiment  is  given 
below  (reference  3-3) . 


Qv 

a 1/2 
1.2(Qv+QL)-0.35(gd)i/^A1 

where  d = inlet-pipe  diameter 

= 0.125  ft.  for  the  plastic  impeller 
= 0.2083  ft.  for  the  bronze  impeller 
and  A^  = inlet-pipe  area. 


(3-2) 


The  details  of  the  void-fraction  calculation  using  the  drift- 
flux  model  are  contained  in  Appendix  D. 

Using  equation  3-1,  for  a water-flow  rate,  QT , an 

Li 

appropriate  air- flow  rate,  Q , was  determined  for  a specified 
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void  fraction.  Thus,  by  setting  appropriate  values  of  Q 

L 

and  Qv,  a void-fraction  range  of  0.0  to  0.62  was  obtained 
during  the  experiment. 

The  acquisition  of  the  two-phase  data  proceeded  in 
much  the  same  manner  as  that  of  the  single-phase  data.  A 
speed  was  set  by  adjusting  the  traction-drive  ratio  of  the 
impeller  motor,  and  checked  by  the  Strobotac.  Using  the  ven- 
turi and  rotameter  calibration  curves,  appropriate  water- 
and  air-flow  rates  were  set  to  yield  desired  void  fractions. 
For  each  value  of  void  fraction  at  different  QT  and  Q the 
static  head  across  the  impeller  was  read  off  the  30-inch 
manometer  and  recorded.  Again,  care  was  taken  to  allow  the 
back  pressure  to  attain  a steady  value,  and  to  remove  any 
air  bubbles  in  the  manometer  lines  before  readings  were 
taken.  Using  the  plastic  impeller  the  above  procedure  was 
repeated  for  190,  250,  318  and  390  rev/min,  and  for  water- 
flow  rates  of  10.0  to  46.0  gal/min,  covering  a void-fraction 
range  of  0.0  to  0.62.  Using  the  bronze  impeller  this  pro- 
cedure was  repeated  for  200,  300,  350  and  400  rev/min,  and  for 
water-flow  rates  of  10.0  to  50.0  gal/min,  covering  a void- 
fraction  range  of  0.0  to  0.40. 

Having  obtained  the  static  head  across  the  pump  from 
the  30- inch-manometer  reading,,  the  dynamic  head  had  to  be 
calculated  so  that  the  total  two-phase  head  could  be  used  in 
the  correlation.  The  calculation  of  the  total  two-phase  head 


46 


was  based  upon  an  effective  two-phase  density  (reference  3-4) , 
p , taking  slip  between  the  vapor  and  liquid  phases  into 
consideration.  The  details  of  the  two-phase  dynamic-head 
calculations  are  given  in  Appendix  C.  The  total  two-phase 
head  across  the  pump  was  then  determined  by  the  following 
equation . 


AH 


otp 


AH.  + 
tp 


, UL2  , 
PL2( (l-a)A2) 


“ Pt  i ( 


Q 


LI 


LI' (l-a)A1< 


2gp 


tp 


where 


tp 


Pvas  + (l-a)pL 
( 1-a)  + as 


Jv(l-a) 

s = 


JLa 


In  this  manner  the  total  head  across  the  pump  was  calculated 
for  63  two-phase  data  points  obtained  from  the  plastic  im- 
peller, and  50  two-phase  data  points  obtained  from  the  bronze 
impeller . 

Having  obtained  a value  of  AHotp  for  each  of  the  two- 
phase  data  points  recorded,  the  head  coefficient,  i|»'  , and 

the  flow  coefficient,  ^tp»  were  calculated  and  recorded  for 
each  data  point.  Lists  of  the  two-phase  data  taken  from  both 
the  plastic  and  bronze  impeller  are  given  in  tables  3-3  and 
3-4  respectively. 
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TABLE  3-3 

TWO-PHASE  PLASTIC- IMPELLER  DATA 


MOTOR 

SPEED 

(rpm) 

TEST 

NO. 

VOID 

FRACTION 

it 

WATER- 
FLOW  RATE 

( ft^/sec) 

AIR- 
FLOW RATE 

( ft'Vsec) 

TOTAL 

HEAD 

(in.  1I20) 

FLOW 

COEFF. 

1' 

HEAD 

COEFF. 

190 

01 

0. 10 

0.0223 

0.0021 

-1.81 

0.0861 

-0.1054 

02 

0.15 

0.0223 

0.0033 

-2.24 

0.0906 

-0.1 305 

03 

0.20 

0.0223 

0.0048 

-2.60 

0.0957 

-0.1514 

04 

0.25 

0.0223 

0.0065 

-4.01 

0.1017 

-0.2335 

05 

0.30 

0.0223 

0.0085 

-4.49 

0.1089 

-0.2615 

06 

0.40 

0.0223 

0.0140 

-5.98 

0.1282 

-0.3483 

07 

0.10 

0.0361 

0.0039 

-3.09 

0.1415 

-0.1700 

08 

0.15 

0.0361 

0.0064 

-4.95 

0.1500 

-0.2883 

09 

0.20 

0.0361 

0.0091 

-6.10 

0.1599 

-0.3553 

10 

0.25 

0.0361 

0.0124 

-6.90 

0.1714 

-0.4019 

250 

18 

0.10 

0.0223 

0.0021 

-1.81 

0.0654 

-0.0609 

19 

0.15 

0.0223 

0.0033 

-2.24 

0.0687 

-0.0754 

20 

0.20 

0.0223 

0 . 0048 

-3.30 

0.0727 

-0.1110 

21 

0.25 

0.0223 

0.0065 

-3.81 

0.0773 

-0.1282 

22 

0.30 

0.0223 

0.0085 

-4.39 

0.0829 

-0.1477 

23 

0.40 

0.0223 

0.0140 

-5.98 

0.0974 

-0.2012 

24 

0.10 

0.0361 

0.0039 

-3.99 

0.1076 

-0.1342 

25 

0.15 

0.0361 

0.0064 

-4.85 

0.1141 

-0.1632 

26 

0.20 

0.0361 

0.0091 

-5.80 

0.1215 

-0.1951 

27 

0.25 

0.0361 

0.0124 

-7.00 

0.1303 

-0.2355 

28 

0.10 

0.0512 

0.0060 

-5.03 

0.1537 

-0.1692 

29 

0.15 

0.0512 

0.0097 

-6.47 

0.1635 

-0.2176 

30 

0.20 

0.0512 

0.0139 

-7.82 

0.1749 

-0.2631 

318 

35 

0.10 

0.0223 

0.0021 

2.19 

0.0515 

0.0455 

36 

0.15 

0.0223 

0.0021 

0.86 

0.0541 

0.0179 

48 


TABLE  3- 3 ( continued ) 


MOTOR 

TEST 

VOID 

WATER- 

AIR- 

TOTAL 

FLOW 

HEAD 

SPEED 

NO. 

FRACTION 

ELOW  RATE 

FLOW  RATE 

HEAD 

COEFF. 

COEFF. 

( « pm) 

< 

( 1 1 ’/see) 

( 1 1 ’/:;.•<) 

(in.  11,0) 

■I1 

<1-* 

3 1 H 

1 7 

0. 20 

0.02  2 ! 

0 . 003  1 

80 

0 . 05  / 2 

-0 . 0582 

38 

0. 25 

0.0223 

0.0065 

-1.91 

0 . 0608 

-0.081  1 

39 

0.30 

0.0223 

0.0D85 

-4 . 59 

0.0651 

-0.0954 

40 

0.40 

0.0223 

0.0140 

-6.18 

0.0766 

-0.1285 

41 

0.10 

0.0361 

0.0039 

0.61 

0.0846 

0.0127 

42 

0.15 

0.0361 

0.0064 

-4.75 

0.0897 

-0.0988 

43 

0.20 

0.0361 

0.0091 

-5.70 

0.0955 

-0.1185 

44 

0.25 

0.0361 

0.0124 

-7.00 

0.1025 

-0.1455 

45 

0.10 

0.0512 

0.0060 

-1.63 

0.1208 

-0.0339 

46 

0.15 

0.0512 

0.0097 

-7.07 

0.1286 

-0.1470 

47 

0.20 

0.0512 

0.0139 

-8.52 

0. 1 375 

-0.1771 

48 

0.  10 

0.0719 

0.0088 

-6.69 

0 . 1 705 

-0. 1 391 

390 

52 

0.10 

0.0223 

0.0021 

2.39 

0.0420 

0.0  3 10 

53 

0.15 

0.0223 

0.0033 

-2.14 

0.0441 

-0.0296 

54 

0.20 

0.0223 

0.0045 

-2.80 

0.0466 

-0.0387 

55 

0.25 

0.0223 

0.0065 

-3.41 

0.0496 

-0.0434 

56 

0.30 

0.0223 

0.0085 

-4.69 

0.0531 

-0.0648 

57 

0.40 

0.0223 

0.014 

-5.88 

0.0625 

-0.0813 

58 

0.10 

0.0361 

0.0039 

2.71 

0.0690 

0.0375 

59 

0.15 

0.0361 

0.0064 

- .05 

0.0731 

-0.0007 

60 

0.20 

0.0361 

0.0091 

-6.20 

0.0779 

-0.0857 

61 

0.25 

0.0361 

0.0124 

-6.90 

0.0835 

-0.0954 

62 

0.10 

0.0512 

0.0060 

. 37 

0.0985 

0.0051 

63 

0.15 

0.0512 

0.0097 

-6.87 

0.1048 

-0.0950 

64 

0.20 

0.0512 

0.0139 

-8.12 

0.1121 

-0.1122 

65 

0.10 

0.0719 

0.0088 

-3.79 

0.1391 

-0.0524 

66 

0.15 

0.0719 

0.0142 

-2.08 

0.1483 

-0. 1670 

4') 

TARLE  1-  3 (com  inu<nl) 


MOTOR 
SPEED 
( rpm) 

TEST 

NO. 

VOID 

FRACTION 

WATER- 
FLOW  RATE 

( f t Vsec) 

AIR- 
FLOW RATE 

( f t Vsec) 

TOTAL 

HEAD 

(in.  II 20 ) 

FLOW 
COEFF . 

4' 

HEAD 

COEFF. 

<l'' 

390 

67 

0. 10 

0.0886 

0.0111 

-4.35 

0.1717 

-0.0602 

190 

69 

0.51 

0.0223 

0.0233 

-8.28 

0.1612 

-0.4822 

250 

70 

0.51 

0.0223 

0.0233 

-8.28 

0.1225 

-0.2785 

71 

0.62 

0.0223 

0.0439 

-14.67 

0.1779. 

-0.4935 

72 

0.37 

0.0361 

0.0233 

-9.52 

0.1596 

-0.3203 

318 

73 

0.51 

0.0223 

0.0233 

-7.78 

0.0963 

-0.1618 

74 

0.62 

0.0223 

0.0439 

-14.27 

0. 1 398 

-0.2967 

75 

0.  37 

0.0361 

0.0233 

-9.52 

0. 1255 

-0.1979 

76 

0.50 

0.0361 

0.0439 

-16.40 

0.1689 

-0.3410 

77 

0.29 

0.0512 

0.0233 

-11.41 

0.1574 

-0.2372 

79 

0.25 

0.0624 

0.0233 

-13.47 

0.1810 

-0.2801 

■30 


tablj:  3-4 


TWO- PHASE  BRONZE- IMPELLER  DATA 


MOTOR 

SPEED 

(rpm) 

TEST 

NO. 

VOID 

FRACTION 

in 

WATER- 
FLOW  RATE 

( ft  Vsec) 

AIR- 
FLOW RATE 

( f t3/sec) 

TOTAL 

HEAD 

(in.  1I20) 

FLOW 

COEFF. 

200 

01 

0.15 

0.0361 

0.0023 

-2.77 

0.1036 

02 

0.20 

0.0361 

0.0033 

-3.16 

0.1063 

03 

0.25 

0.0361 

0.0044 

-3.86 

0.1092 

04 

0.10 

0.0512 

0.0035 

-4.54 

0.1475 

05 

0.15 

0.0512 

0.0056 

-10.43 

0.1532 

06 

0.20 

0.0512 

0.0080 

-12.42 

0.1597 

300 

07 

0.15 

0.0361 

0.0023 

-0.97 

0.0691 

08 

0.20 

0.0361 

0.0033 

-1.96 

0.0709 

09 

0.25 

0.0361 

0.0044 

-2.66 

0.0728 

10 

0.10 

0.0512 

0.0035 

-3.84 

0.0984 

11 

0.15 

0.0512 

0.0056 

-10.73 

0.1021 

12 

0.20 

0.0512 

0.0080 

-11.32 

0.1065 

13 

0.10 

0.0719 

0.0063 

-12.08 

0.1406 

14 

0.15 

0.0719 

0.0101 

-23.86 

0.1475 

15 

0.10 

0.0886 

0.0086 

-16.52 

0.1748 

350 

16 

0.15 

0.0361 

0.0023 

-0.47 

0.0592 

17 

0.20 

0.0361 

0.0033 

-0.76 

0.0607 

18 

0.25 

0.0361 

0.0044 

-0.96 

0.0624 

19 

0.10 

0.0512 

0.0035 

-3.44 

0.0843 

20 

0.15 

0.0512 

0.0056 

-9.93 

0.0875 

21 

0.20 

0.0512 

0.0060 

-10.92 

0.0912 

22 

0.  10 

0.0719 

0.0063 

-8.28 

0.1205 

23 

0.15 

0.0719 

0.0101 

-23.36 

0.1264 

24 

0.10 

0.0886 

0.0086 

-10.82 

0.1498 

25 

0.  10 

0. 1022 

0.0104 

-21.55 

0.1736 

HEAD 
COEFF . 

V 

-0.1825 

-0.2082 

-0.2544 

-0.2992 

-0.6873 

-0.8184 

-0.0284 

-0.0574 

-0.0779 

-0.1125 

-0.3142 

-0.3315 

-0.3538 

-0.6988 

-0.4838 

-0.0101 

-0.0164 

-0.0207 

-0.0740 

-0.2137 

-0.2350 

-0.1782 

-0.5026 

-0.2328 

-0.4637 
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TABLE  1-4 (continued ) 


MOTOR 
SPEED 
( rpni) 

TEST 

NO. 

VOID 

FRACTION 

WATER- 
FLOW  RATE 

( ft  Vs  PC) 

AIR- 
FLOW RATE 

( f t J/spc) 

TOTAL 

HEAD 

(in.  1 1 ^ 0 ) 

FLOW 
COEFF . 

'!• 

HEAD 

COEFF. 

400 

26 

0.15 

0.0361 

0.00  2 1 

0.2  1 

0.0518 

0.00  18 

27 

0.20 

0.0361 

0.003  1 

-0. 26 

0.05  1 1 

-0.004  1 

28 

0.25 

0.0361 

0 . 0044 

-0.36 

0.0546 

-0.0059 

29 

0.10 

0.0512 

0.0035 

-3.24 

0.0738 

-0.0534 

30 

0.15 

0.0512 

0.0056 

-7.23 

0.0766 

-0.1191 

31 

0.20 

0.0512 

0.0080 

-7.72 

0.0798 

-0.1272 

32 

0.10 

0.0719 

0.0063 

-5.08 

0.1055 

-0.0837 

33 

0.15 

0.0719 

0.0101 

-21.66 

0.1106 

-0. 3568 

34 

0.10 

0.0886 

0.0086 

-5.82 

0.1311 

-0.0959 

35 

0.10 

0.1022 

0.0104 

-15.55 

0.1519 

-0.2562 

200 

36 

0.30 

0.0361 

0.0058 

-4.15 

0.1130 

-0.2735 

37 

0.40 

0.0  16 1 

0.0095 

-4 . 72 

0. 1 2 10 

-0.1110 

38 

0. 10 

0.0624 

0.0050 

-8.11 

0.  1 8 1 H 

-0.5  14  4 

300 

40 

0. 30 

0.0361 

0.0058 

-2.85 

0.0753 

-0.08  15 

41 

0.40 

0.0361 

0.0095 

-4.02 

0.0820 

-0.1177 

42 

0.10 

0.0624 

0.0050 

-5.81 

0.1212 

-0.1702 

43 

0.20 

0.0624 

0.0116 

-19.98 

0.1331 

-0.5851 

350 

44 

0.30 

0.0361 

0.0058 

-1.25 

0.0646 

-0.0269 

4 5 

0.40 

0.0361 

0.0095 

-1.92 

0.0703 

-0.0413 

46 

0.10 

0.0624 

0.0050 

-4.71 

0.1039 

-0.1013 

47 

0.20 

0.0624 

0.0116 

-17.48 

0.1141 

-0. 3761 

400 

48 

0.  30 

0.0361 

0.0058 

-0.  35 

0.0565 

-0.0058 

49 

0.40 

0.0361 

0.0095 

-1.02 

0.0615 

-0.0168 

50 

0.10 

0.0624 

0.0050 

-3.51 

0.0909 

-0.0578 

51 

0.20 

0.0624 

0.0116 

-14.58 

0.0998 

-0.2402 

52 


CHAPTER  4 

DERIVATION  OF  THEORETICAL  CHARACTERISTICS 

SINGLE  PHASE 

The  first  step  in  determining  the  theoretical  single- 
phase pump  characteristic  was  to  apply  Euler's  turbomachinery 
equation  to  express  the  change  in  enthalpy  of,  or  work  done 
on,  the  fluid  flowing  through  the  pump,  to  the  change  in 
moment  of  momentum  of  the  rotor  inlet  and  outlet  flows.  In 
applying  Euler's  equation,  one-dimensional  steady  flow  was 
assumed.  Euler's  equation  for  adiabatic  flow  can  be  expressed 
mathematically  as  follows. 

- 9c<ho2  - hol>  * = U2C92  - U1CQ1  (4-1) 

Since  a f low-straightener  was  used  in  the  inlet  piping  and 
inlet  swirl  vanes  were  not  used  throughout  the  experiment, 

C0^  = 0 and  equation  4-1  can  be  simplified  to 

g Ah  = u„CQO 

^c  o 2 02 

For  single-phase  flow,  the  simple  (one-dimensional) 
velocity  triangles  at  the  rotor  tips  of  the  plastic  and 
bronze  impellers  are  shown  in  figures  4-1  and  4-2  respective- 
ly. The  angle  is  the  relative  flow  angle  at  the  rotor- 
outlet  diameter,  and  is  equax  to  the  rotor  blade  angle, 
less  the  slip  angle,  6.  The  slip  angle  for  each  of  the  im- 
pellers was  calculated  using  the  Busemann  slip-factor 
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a 

2 


FIGURE  4-1.  VELOCITY  DIAGRAM  AT  TIIE  PLASTIC  IMPELLER'S 
BLADE  TIP 
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correlation.  The  details  of  these  calculations  can  be 
found  in  Appendix  A.  The  results  of  these  calculations  were 
B2  = 41.76°  for  the  plastic  impeller,  and  $2  = 12.5°  for  the 
bronze  impeller. 

From  the  geometry  of  the  velocity  triangles 
Cm2 


C92  U2  ” tan  £. 


Therefore, 


g Ah  = u." 
yc  o 2 


U2Cm2 
tan  B. 


(4-2) 


Dividing  through  equation  4-2  by  u2  yields 

C 


g Ah 
^c  o 


u. 


= 1.0  - 


'm2 


u. 


tan  B. 


Now,  for  the  specific  case  of  loss-less  incompressible 
flow,  the  theoretical  total-enthalpy  change  is  proportional 
to  the  theoretical  total-head  change.  Thus,  ths  expression 
for  the  theoretical  single-phase  pump  characteristic  is 


gAH  g Ah  C _ 

^ o _ ~c  o __  i a m2  .. 

u 2 u 2 U2 


tan  B. 


(4-3) 


2 2 

Now,  employing  the  definitions  of  flow  coefficient, 
C 

* = -i" 

T u 

and  head  coefficient, 
gAH 

. i - o 

r = — 

U„ 
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equation  4-3  becomes 

^2 

^ spth  tan  fl2 


(4-4) 


For  the  plastic  impeller,  using  Busemann's  slip-factor 
correlation  and  the  geometry  of  the  impeller,  the  theoretical 
single-phase  characteristic  curve,  in  the  form  of  equation 
4-4,  was  determined  to  be 

tp  * spth  = 1,0  " 1-120014  <J>2  (4-5) 

Following  the  same  procedure,  the  theoretical  single-phase 
characteristic  curve  equation  for  the  bronze  impeller  was 
calculated  to  be 

♦'spth-  1-°  - 4-511  *2  l4-6> 

TWO  PHASE 

The  Euler  equation  is  applicable  to  both  single-  and 
two-phase  flow.  However,  in  two-phase  flow  there  are  two 
separate  mass  flows  with  presumably  different  tangential 
and  radial  velocities  at  both  rotor  inlet  and  outlet.  Euler's 
equation  for  two-phase  flow  with  zero  inlet  swirl  then 
becomes 


*'tpth  - (1-x2) 


'6L2 

u„ 


+ X. 


'6v2 

u„ 


(4-7) 


From  the  geometry  of  the  velocity  triangles  at  the  rotor 
outlet,  equation  4-7  can  be  transformed  to 
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tpth 


(l-x2) [1.0 


^mL2 

(u2tan  e2) 


] + x2[1.0  - 


Cmv2 

(u2tan  32) 


where  the  vapor  and  liquid  phases  are  assumed  to  leave  the 
rotor  at  the  same  relative  angle.  There  may  in  fact  be  some 
difference  in  their  direction,  but  the  effect  will  be  small 
(reference  4-1) . 

Substituting  for 


'mL2 


lhL2/(pL2AL2) 


*Td-x2)/ 


^PL2  ^1-a2^ AT2^ 


Cmv2  lflv2//(pv2Av2)  lflTX2/ [ pv2a2AT2 J 

and  defining  a two-phase  flow  coefficient 

»tP  s ss*  “ V<ptpAu) 

where  the  mean  two-phase  density  is  defined  as 
ptp  5 (1"a)pL  + apv 

then 


= 1.0 


- 2**-  [i.o  + tA  ^1 


tpth  J'*'/  tan  ' (l-o^)  pl2 


t(l-x  »2  + —3-2-  22  x 2] 
z a2  pv2  1 


Defining  a two-phase  flow  function 
a = (a/  (1-a) 
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a slip  velocity  ratio, 

s = C /CT 
v L 


and  quality, 

x = as/(l  + as) 

it  can  be  shown  (reference  4-2)  that  the  theoretical  two- 
phase  characteristic  equation  for  a centrifugal  pump  with 
zero  inlet  swirl  becomes 


tpth 


= 1.0  - f 


Jtp2 


tp2  tan  8 2 


(4-8) 


where 


_ (1  + a)  (1  + as  ) 


tp 


(1  4-  as) 

Equation  4-8  is  applicable  to  both  single-  and  two- 
phase  flow.  In  single-phase  flow  a = 0 = a,  and  ffc^  = 1.0, 
transforming  equation  4-8  back  into  equation  4-6  previously 
developed  for  single-phase  flow. 

Substituting  the  appropriate  values  for  all  of  the  two- 
phase  data  points  taken,  using  each  of  the  impellers,  into  a 
computer  program  produced  by  Tak-Chee  Stephen  Chan  (Appen- 
dix E)  yielded  the  theoretical  two-phase  head  coefficients 
for  each  of  the  selected  data  points. 

Thus,  having  calculated  the  theoretical  single-  and  two- 
phase  characteristic  curve  equations  for  both  the  plastic  and 
bronze  impellers,  the  next  step  was  to  determine  the  head-loss 


ratio. 
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CHAPTER  5 

CALCULATION  OF  THE  HEAD-LOSS  RATIO 


The  head- loss  ratio,  H*,  the  ratio  of  the  flow  losses 
in  two-phase  flow  to  those  in  single-phase  flow,  was  calcu- 
lated so  that  it  could  be  plotted  versus  void  fraction  and 
the  results  compared  to  the  correlation  proposed  by  Mikiel- 
ewicz  and  Wilson.  Mathematically,  the  head-loss  ratio  can 
be  defined  as 


H*  = !_t£th ^_t£ 

■ *’spth  - *'sp 


(5-1) 


H*  for  every  two-phase  data  point  was  then  determined 
in  the  following  manner.  Values  of  ' fcp  were  determined  from 
the  definition  of  the  head  coefficient,  the  two-phase  total- 
head  data,  and  the  impeller  geometry.  For  the  plastic  im- 
peller, values  of  i}*'  were  obtained  from  figure  3-1  for 


appropriate  values  of 


For  the  bronze  impeller,  values  of 


[f)’  were  calculated  by  substituting  appropriate  values  of  <f> 
sp 

into  equation  3-1.  Finally,  values  of  and  ^ ' spth 

were  calculated  using  equations  4-8  and  4-4  respectively,  for 
appropriate  values  of  flow  coefficient  and  82* 

The  values  of  head-loss  ratio,  void  fraction,  and  the 


associated  head  coefficients  used  in  determining  H*,  for  the 
data  points  taken  using  the  plastic  impeller,  are  given  in 
table  5-1.  Table  5-2  contains  similar  results  from  the  bronze- 
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impeller  tests.  The  values  contained  in  table  5-2  were 
arrived  at  through  the  use  of  a computer  program  developed 
by  Tak-Chee  Stephen  Chan  (Appendix  E) . 
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TABLE  5-1 


CALCULATION  OF  HEAD-LOSS  RATIO 
FOR  THE  PLASTIC  IMPELLER  DATA 


TEST 

NO. 

a 

<P 

^ spth 

1 1>' 

sp 

*'tpth 

iK . 
tp 

H* 

01 

0.10 

0.0861 

0.9036 

0.2150 

0.9036 

-0.1054 

1.4653 

02 

0.15 

0.0906 

0.8985 

0.2070 

0.8986 

-0.1305 

1.4881 

03 

0.20 

0.0957 

0.8928 

0.1988 

0.8928 

-0.1514 

1.5046 

04 

0.25 

0.1017 

0.8861 

0.1880 

0.8860 

-0.2335 

1.6037 

05 

0.30 

0.1089 

0.8780 

0.1750 

0.8780 

-0.2615 

1.6208 

06 

0.40 

0.1282 

0.8564 

0.1400 

0.8564 

-0.3483 

1.6815 

07 

0.10 

C 1415 

0.8415 

0.1135 

0.8415 

-0.1700 

1.4031 

08 

0.15 

0.1500 

0.8320 

0.0945 

0.8319 

-0.2883 

1.5189 

09 

0.20 

0.1599 

0.8209 

0.0724 

0.8209 

-0.3553 

1.5714 

10 

0.25 

0.1714 

0.8080 

0.0375 

0.8080 

-0.4019 

1.5702 

10 

0.10 

0.0654 

0.9268 

0.2465 

0.9267 

-0.0609 

1.4518 

19 

0.15 

0 . 0687 

0.9231 

0.2413 

0.9229 

-0.0754 

1.4643 

20 

0.20 

0.0727 

0.9186 

0.2355 

0.9185 

-0.1110 

1.5072 

21 

0.25 

0.0773 

0.9134 

0.2290 

0.9134 

-0.1282 

1.5219 

22 

0.25 

0.0773 

0.9134 

0.2290 

0.9134 

-0.1282 

1.5219 

23 

0.30 

0.0828 

0.9073 

0.2195 

0.9073 

-0.1477 

1.5339 

24 

0.10 

0.1076 

0.8795 

0.1780 

0.8795 

-0.1342 

1.4451 

25 

0.15 

0.1141 

0.8722 

0.1715 

0.8723 

-0.1632 

1.4777 

26 

0.20 

0.1215 

0.8639 

0.1533 

0.8639 

-0.1951 

1.4902 

27  0.25  0.1303  0.8541  0.1375  0.8540  -0.2355  1.5205 
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TABLE  5-1  (continued) 


TEST 

NO. 

a 

4> 

^ spth 

ip ' 

sp 

*’tpth 

ip ' 

W tp 

28 

0.10 

0.1537 

0.8279 

0.0850 

0.8279 

-0.1692 

29 

0.  15 

0.1635 

0.8169 

0.0625 

0.8168 

-0.2176 

30 

0.20 

0.1749 

0.8041 

0.02*5 

0.8041 

-0.2631 

35 

0.10 

0.0515 

0.9423 

0.2675 

0.9424 

0.0455 

36 

0.15 

0.0541 

0.9394 

0.2635 

0.9394 

0.0179 

37 

0.20 

0.0572 

0.9359 

0.2590 

0.9360 

-0.0582 

38 

0.25 

0.0608 

0.9319 

0.2535 

0.9319 

-0.0813 

39 

0.30 

0.0651 

0.9271 

0.2470 

0.9271 

-0.0954 

40 

0.40 

0.0766 

0.9142 

0.2310 

0.9142 

-0.1285 

41 

0.10 

0.0846 

0.9053 

0.2170 

0.9053 

0.0127 

42 

0.15 

0.0897 

0.8995 

0.2085 

0.8996 

-0.0988 

43 

0.20 

0.0955 

0.8930 

0.1990 

0.8930 

-0.1185 

44 

0.25 

0.1025 

0.8852 

0.1870 

0.8853 

-0.1455 

45 

0.10 

0.1208 

0.8647 

0.1565 

0.8647 

-0.0339 

46 

0.15 

0.1286 

0.8560 

0.1400 

0.8560 

-0.1470 

47 

0.20 

0.1375 

0.8460 

0.1220 

0.8460 

-0.1771 

48 

0.10 

0.1705 

0.8090 

0.0425 

0.8090 

-0.1391 

52 

0.10 

0.0420 

0.9530 

0.2825 

0.9530 

0.0330 

53 

0.15 

0.0441 

0.9506 

0.2800 

0.9506 

-0.0296 

54 

0.20 

0.0466 

0.9478 

0.2770 

0.9478 

-0.0387 

55 

0.25 

0.0496 

0.9445 

0.2710 

0.9445 

-0.0434 

H* 

1.3422 

1.3712 

1.3689 

1.3291 

1.3633 

1.4686 

1.4935 

1.5035 

1.5262 

1.2969 

1.4447 

1.4574 

1.4763 

1.2688 

1.4009 

1.4131 

1.2369 

1.3722 

1.4616 

1.4706 

1.4669 


56 


0.30  0.0531  0.9405  0.2650  0.9406  -0.0648  1.4883 
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TABLE  5-1  (continued) 


TEST 

NO. 

a 

^ spth 

• 

sp 

*'tpth 

ip ' 

V tp 

57 

0.40 

0.0625 

0.9300 

0.2510 

0.9301 

-0.0813 

59 

0.15 

0.0731 

0.9181 

0.2350 

0.9181 

-0.0007 

60 

0.20 

0.0779 

0.9128 

0.2275 

0.9128 

-0.0857 

61 

0.25 

0.0835 

0.9065 

0.2190 

0.9064 

-0.0954 

62 

0.10 

0.0985 

0.8897 

0.1940 

0.8897 

0.0051 

63 

0.15 

0.1048 

0.8826 

0.1830 

0.8826 

-0.0950 

64 

0.20 

0.1121 

0.8745 

0.1710 

0.8744 

-0.1122 

65 

0.10 

0.1391 

0.8442 

0.1185 

0.8443 

-0.0524 

66 

0.15 

0.1483 

0.8339 

0.0975 

0.8339 

-0.1670 

67 

0.10 

0.1717 

0.8077 

0.0375 

0.8077 

-0.0602 

69 

0.51 

0.1612 

0.8195 

0.0690 

0.8234 

-0.4822 

70 

0.51 

0.1225 

0.8628 

0.1515 

0.8628 

-0.2785 

71 

0.62 

0.1779 

0.8008 

0.0075 

0.8007 

-0.4935 

72 

0.37 

0.1596 

0.8213 

0.0725 

0.8212 

-0.3203 

73 

0.51 

0.0963 

0.8921 

0.1975 

0.8921 

-0.1618 

74 

0.62 

0.1398 

0.8434 

0.1175 

0.8434 

-0.2967 

75 

0.37 

0.1255 

0.8594 

0.1463 

0.8595 

-0.1979 

76 

0.50 

0.1689 

0.8108 

0.0470 

0.8107 

-0.3410 

H* 

1.4895 

1.3450 

1.4571 

1.4573 

1.2715 

1.3973 

1.4025 

1.2356 

1.3592 

1.1268 

1.7400 

1.6045 

1.6316 

1.5246 

1.5172 

1.5705 

1.4827 

1.5078 


77  0.29  0.1574  0.8237  0.0780  0.8237  -0.2372  1.4227 
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TABLE  5-2 


CALCULATION  OF  HEAD- LOSS  RATIO 
FOR  THE  BRONZE- IMPELLER  DATA 


TEST 

NO. 

a 

<f> 

^ spth 

ip ' 

sp 

^ tpth 

ib' 

V tp 

H* 

01 

0.15 

0.1036 

0.5328 

0.1950 

0.5328 

-0.1825 

2.1173 

02 

0.20 

0.1063 

0.5206 

0.1800 

0.5206 

-0.2082 

2.1400 

03 

0.25 

0.1092 

0.5072 

0.1629 

0.5072 

-0.2544 

2.2118 

04 

0.10 

0.1475 

0.3345 

-0.1324 

0.3345 

-0.2992 

1.3571 

05 

0.15 

0.1532 

0.3089 

-0.1876 

0.3089 

-0.6873 

2.0065 

06 

0.20 

0.1597 

0.2797 

-0.2532 

0.2797 

-0.8184 

2.0607 

07 

0.15 

0.0691 

0.6885 

0.3391 

0.6885 

-0.0284 

2.0520 

08 

0.20 

0.0709 

0.6804 

0.3335 

0.6804 

-0.0574 

2.1270 

09 

0.25 

0.0728 

0.6715 

0.3271 

0.6715 

-0.0779 

2.1762 

10 

0.10 

0.0984 

0.5563 

0.2221 

0.5563 

-0.1125 

2.0012 

11 

0.15 

0.1021 

0.5393 

0.2027 

0.5393 

-0.3142 

2.5356 

12 

0.20 

0.1065 

0.5198 

0.1790 

0.5198 

-0.3315 

2.4981 

13 

0.10 

0.1406 

0.3657 

-0.0686 

0.3657 

-0.3538 

1.6565 

14 

0.15 

0.1475 

0.3349 

-0.1316 

0.3349 

-0.6988 

2.2160 

15 

0.10 

0.1748 

0.2116 

-0.4102 

0.2116 

-0.4838 

1.1183 

16 

0.15 

0.0592 

0.7330 

0.3666 

0.7330 

-0.0101 

2.0279 

17 

0.20 

0.0607 

0.7261 

0.3627 

0.7261 

-0.0164 

2.0429 

18 

0.25 

0.0624 

0.7184 

0. 3582 

0.7184 

-0.0207 

2.0517 

19 

0.10 

0.0843 

0.6197 

0.2853 

0.6197 

-0.0740 

2.0743 

20  0.15  0.0875  0.6051  0.2719  0.6051  -0.2137  2.4576 
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TABLE  5-2  (continued) 


TEST 

NO. 

a 

4> 

^ spth 

i|>' 

sp 

* ' tpth 

ib  ' 

V tp 

H* 

21 

0.20 

0.0912 

0.5884 

0.2558 

0.5884 

-0.2350 

2.4758 

22 

0.10 

0.1205 

0.4563 

0.0905 

0.4563 

-0.1782 

7346 

23 

0.15 

0.1264 

0.4299 

0.0482 

0.4299 

-0.5026 

2.4434 

24 

0.10 

0.1498 

0.3242 

-0.1543 

0.3242 

-0.2328 

1.1641 

25 

0.10 

0.1736 

0.2171 

-0.3976 

0.2171 

-0.4637 

1.1076 

26 

0.15 

0.0518 

0.7664 

0.3834 

0.7664 

0.0038 

1.9912 

27 

0.20 

0.0531 

0.7603 

0.3806 

0.7603 

-0.0043 

2.0135 

28 

0.25 

0.0546 

0.7536 

0.3773 

0.7536 

-0.0059 

2.0186 

29 

0.10 

0.0738 

0.6672 

0.3240 

0.6672 

-0.0534 

2.0994 

30 

0.15 

0.0766 

0.6545 

0.3143 

0.6545 

-0.1191 

2.2739 

31 

0.20 

0.0798 

0.6399 

0.3025 

0.6399 

-0.1272 

2.2739 

32 

0.10 

0.1055 

0.5243 

0.1846 

0.5243 

-0.0837 

1.7897 

33 

0.15 

0.1106 

0.5012 

0.1549 

0.5012 

-0.3568 

2.4777 

34 

0.10 

0.1311 

0.4087 

0.0118 

0.4087 

-0.0959 

1.2714 

35 

0.10 

0.1519 

0.3150 

-0.1742 

0.3150 

-0.2562 

1.1675 

36 

0.30 

0.1130 

0.4902 

0.1400 

0.4902 

-0.2735 

2.1807 

37 

0.40 

0.1230 

0.4452 

0.0731 

0.4452 

-0.3110 

2.0326 

38 

0.10 

0.1818 

0.1799 

-0.4800 

0.1799 

-0.5344 

1.0824 

40 

0.30 

0.0753 

0.6601 

0.3187 

0.6601 

-0.0835 

2.1775 

41 

0.40 

0.0820 

0.6301 

0.2944 

0.6301 

-0.1177 

2.2273 

42 

0.10 

0.1212 

0.4533 

0.0859 

0.4533 

-0.1702 

1.6969 

43  0.20  0.1331  0.3998  -0.0041  0.3998  -0.5851  2.4385 
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TABLE  5-2  (continued) 


TEST 

NO. 

a 

0/ ' 

v spth 

i|»' 

sp 

^ tpth 

r 

H* 

44 

0.30 

0.0646 

0.7087 

0.3523 

0.7087 

-0.0269 

2.0638 

45 

0.40 

0.0703 

0.6830 

0.3353 

0.6830 

-0.0413 

2.0833 

46 

0.10 

0.1039 

0.5314 

0.1933 

0.5314 

-0.1013 

1.8713 

47 

0.20 

0.1141 

0.4855 

0.1335 

0.4855 

-0.3761 

2.4475 

48 

0.30 

0.0565 

0.7451 

0.3730 

0.7451 

-0.0058 

2.0181 

49 

0.40 

0.0615 

0.7226 

0.3607 

0.7226 

-0.0168 

2.0428 

50 

0.10 

0.0909 

0.5900 

0.2574 

0.5900 

-0.0578 

1.9478 

51 

0.20 

0.0998 

0.5498 

0.2148 

0.5498 

-0.2402 

2.3583 

67 


CHAPTER  6 

RESULTS 

HEAD-LOSS  RATIO  VERSUS  VOID  FRACTION 

Having  calculated  H*  for  all  the  data  points  taken,  the 
head- loss  ratio  was  then  plotted  versus  void  fraction. 

Figure  6-1  is  a plot  of  H*  versus  a with  flow  coefficient  as 
a parameter,  for  the  data  taken  using  the  plastic  impeller. 

No  sharp  rise  in  the  head- loss  ratio  was  observed  at  0.20 
void  fraction,,  as  was  the  case  for  the  data  from  the  Olson 
experiments  (reference  3-4) . Rather,  a gradual  increase  in 
H*  with  a was  shown  with  some  correlation  of  increasing  H* 
with  increasing  flow  coefficient  at  void  fraction  of  0.25, 
and  flow  coefficients  of  <_  0.13.  Below  a void  fraction  of 
0.25  there  was  too  much  scatter  in  the  data  to  notice  any  cor- 
relation of  H*  with  4>tp»  while  at  flow  coefficient  of  0.13 
H*  tended  to  decrease  with  increasing  <(>  . The  fact  that  the 

maximum  value  achieved  for  H*  was  less  than  1.8  indicated  that 
the  plastic  impeller  was  very  inefficient  in  single-phase 
operation,  and  the  introduction  of  two-phase  flow  did  not  make 
the  performance  deteriorate  as  greatly  as  it  would  have  if  the 
impeller  had  had  a high  single-phase  efficiency. 

Figure  6-2  shows  the  plots  of  *'tpth,  ^'spth'  ^'tp 
and  '!>’Sp  versus  flow  coefficient  for  the  data  and  geometry  of 
the  bronze  impeller.  Calculating  H*,  for  appropriate  values 


HEAD-LOSS  RATIO  , 


1.5 

1.4 


VOID  FRACTION,  d 

FIGURE  6-1.  HEAD- LOGS  RATIO  VS.  VOID  FRACTION  FOR  PLASTIC- 
IMPELLER  DATA  WITH  FLOW  COEFFICIENT  AS  A 
PARAMETER 
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FIGURE  6-2.  ACTUAL  AND  THEORETICAL  CHARACTERISTICS  OF 

SINGLE-  AND  TWO-PHASE  TESTS  FOR  BRONZE- IMPELLER 
DATA 
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of  flow  coefficient,  from  these  curves  and  plotting  versus 
void  fraction  yielded  figure  6-3.  Figure  6-3  again  shows  a 
gradual  increase  of  H*  with  a with  no  sudden  jump  at 
a = 0.20.  To  better  determine  the  existence  of  any  possible 

I 

correlation  of  H*  with  flow  coefficients,  plots  of  H*  versus 

a for  different  ranges  of  flow  coefficient  were  made  and 

are  given  in  figures  6-4  through  6-7.  These  plots  show  a 

better  correlation  of  increasing  H*  with  increasing  flow 

coefficient  than  was  evident  from  the  plastic-impeller  data 

plots.  For  the  bronze  impeller  H*  increased  with  increasing 

dK  for  void  fractions  of  > 0.15,  and  for  flow  coefficients 
tp 

of  £ 0.123.  Above  4>tp  = 0.123  the  head-loss  ratio  again 

decreased  as  tp.  continued  to  increase.  The  maximum  value 
tp 

achieved  by  H*  for  the  bronze-impeller  data  was  2.6,  indicat- 
ing that  the  bronze  impeller  is  more  efficient  in  single-phase 
operation  than  is  the  plastic  impeller,  as  was  expected. 

FLOW  VISUALIZATION 

Throughout  the  testing  careful  observations  of  the  flow 
patterns  in  the  inlet  pipe,  the  impeller  passages,  and  the 
scroll  were  made.  Using  the  plastic  impeller  the  major  sig- 
nificant observation  was  the  growth  of  a large  air  cavity  in 
the  impeller  eye  and  its  expansion  across  the  blades  into  the 
volute  as  the  void  fraction  was  increased.  At  void  fractions 
of  0.10  to  0.15  the  impeller  was  covered  and  filled  with 
small  air  bubbles  recirculating  through  the  impeller.  As  the 
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VOID  FRACTION 


FIGURE  6-4.  HEAD-LOSS  RATIO  VS.  VOID  FRACTION  FOR  BRONZE- 
IMPELLER  DATA  WITH  0.0  < <fr  < 0.076 
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VOID  FRACTION 

FIGURE  6-5.  IIEAD-LOSS  RATIO  VS.  VOID  FRACTION  FOR  BRONZE- 
IMPELLER  DATA  WITH  0.076  < <)>  < 0.102 
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FIGURE  6-6.  IIEAD-LOGG  RATIO  VG.  VOID  FRACTION  FOR  BRONZE- 
IMPELLER  DATA  WITH  0.102  < <j>  < 0.140 
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void  fraction  was  increased  further,  this  air  cavity  in  the 
impeller  eye  grew  in  size  across  the  blades  until  at  approx- 
imately a = 0.20  the  cavity  had  expanded  out  to  the  blade 
tips,  with  water  flowing  only  along  the  bottom  of  the  im- 
peller. As  the  air  cavity  expanded  out  to  the  blade  tips  a 
corresponding  sudden  worsening  of  the  head  degradation  was 
noticed  on  the  static-head  manometer,  connected  across  the 
pump.  As  the  void  fraction  increased  even  further,  the  air 
cavity  grew  into  the  volute  with  the  head  rise  continuing 
to  decrease,  but  degrading  at  a slower  rate  than  that  exper- 
ienced while  the  air  cavity  grew  to  the  blade  tips.  Once 
this  large  air  cavity  was  formed,  the  inlet  void  fraction 
could  be  reduced  to  zero  (single-phase  water)  without  immed- 
iately affecting  the  cavity  size.  It  took  from  10  to  20 
seconds  for  the  cavity  to  collapse  and  for  the  pump  to  return 
to  "normal"  single-phase  operation.  In  addition,  this  air 
cavity  exhibited  an  unsteady  bi-stable  nature  at  void  fract- 
ions between  0.17  and  0.20  by  successively  popping  out  to  the 
blade  tips  then  collapsing  back  into  the  impeller  with  a 
random  frequency.  Photos  6-1  through  6-4  and  6-5  through  6-8 
show  the  growth  of  this  air  cavity  from  a void  fraction  of 
0.10  to  0.40  and  from  0.10  to  0.25  respectively. 

The  flow  in  the  inlet  pipe  was  also  observed  during 
the  plastic-impeller  tests.  Determination  of  the  void  fract- 
ion from  visual  observations  was  found  to  be  virtually 


PHOTO  6-3.  PLASTIC  IMPELLER  DURING  OPERATION  AT  190  RPM 
a - 0.25  AND  <})  - 0.1715 
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PLASTIC  IMPELLER  DURING  OPERATION  AT  250  RPM 
a « 0.40  AND  <J>  = 0.0593 


PHOTO  6-5.  PLASTIC  IMPELLER  DURINC  OPERATION  AT  31G  RPM 
a = 0.10  AND  <J>  ~ 0.06  39 


PHOTO  6-6.  PLASTIC  IMPELLER  DURING  OPERATION  AT  390  RPI1 
a = 0.15  NAD  <p  - 0.0731 


00 


PHOTO 


PHOTO 


-8.  PLASTIC  IMPELLER  DURING  OPERATION  AT  390  RPM, 
a = 0.25  AND  $ - 0.0035 
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impossible,  although  a pulsing  slug-flow  was  seen  to  develop 
at  a void  fraction  of  0.20  - 0.25. 

Since  the  bronze  impeller  was  shrouded  no  observations 
of  the  flow  within  the  blade  channels  was  possible.  However, 
the  flow  in  the  inlet  pipe,  at  the  blade  tips,  and  in  the 
volute  was  carefully  noted.  Photos  6-9  through  6-12  show  the 
inlet-pipe  flow  during  void  fractions  of  0.10  to  0.40.  For 
void  fractions  of  0.10  - 0.20  the  inlet  flow  was  filled  with 
dispersed  buboles  traveling  down  the  pipe  with  the  water  flow. 
However,  at  void  fraction  of  0.20  the  bubbles  became  more 
mixed  and,  in  various  sections  of  the  pipe,  reversed  their  down- 
ward flow  direction.  As  the  void  fraction  was  increased  from 
0.20  to  0.40  this  flow  reversal  became  stronger  and  a pulsing 
slug-flow  developed. 

No  expanding  air  cavity  was  noticed  at  the  blade  tips 
or  in  the  scroll  during  the  bronze-impeller  tests.  Through- 
out the  range  of  void  fraction  the  impeller  periphery  was  sur- 
rounded with  small  air  bubbles.  The  size  of  these  bubbles  in- 
creased with  a decrease  in  the  pump  rotational  speed  and  as 
the  bubbles  flowed  from  the  narrower  to  the  wider  scroll  area. 
In  addition,  a stagnant  water  and  air  void  filled  from  one- 
quarter  to  one-half  of  the  scroll  exit  area  adjacent  to  the 
cutwater  throughout  the  two-phase  testing.  Typical  flow  con- 
figurations at  the  blade  tips  and  in  the  scroll  are  shown  in 
photos  6-13  through  6-20. 


PHOTO  6-15.  BRONZE  IMPELLER  DURING  OPERATION  AT  400  RPM, 
a = 0.10  AND  (J>  =--  0.0738 


PHOTO  6-16.  BRONZE  IMPELLER  DURING  OPERATION  AT  400  RPM, 
a = 0.20  AND  <J>  - 0.0/99 
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PHOTO  6-17.  BRONZE  IMPELLER  DURING  OPERATION  AT  200  RPM, 
a = 0.30  AND  <J>  = 0.1131 


PHOTO  6-18.  BRONZE  IMPELLER  DURING  OPERATION  AT  200  RPM, 
a = 0.40  AND  <J>  = 0.1230 


PHOTO  6-19.  BRONZE  IMPELLER  DURING  OPERATION  AT  400  RPM, 
a =■■  0.30  AND  <J>  = 0.0565 


PHOTO  6-20.  BRONZE  IMPELLER  DURING  OPERATION  AT  400  RPM, 
a = 0.40  AND  <p  = 0.0615 
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CHAPTER  7 

CONCLUSIONS 

The  most  significant  conclusion  drawn  from  these  ex- 
periments is  the  correlation  of  H*  with  flow  coefficient  as 
well  as  with  void  fraction.  Although  this  correlation  is 
visible  only  for  a >_  0.25  in  the  plastic-impeller  data,  the 
data  for  the  bronze  impeller  show  a much  stronger  correlation 
for  a >_  0.15.  The  lack  of  any  definite  correlation  between 
H*  and  flow  coefficient  at  the  lower  void  fractions  may  be 
due  to  two  possibilities.  First  of  all,  the  flow  regime  of 
the  inlet  flow  for  0.10  _<  a <_  0.20  may  be  continually  chang- 
ing from  dispersed  bubble  to  elongated  bubble  to  slug-flow 
and  back  again  with  no  apparent  regularity.  Also,  for  low 
void  fractions  at  flow  coefficients  below  design  point,  H* 
increases  with  increasing  <f>  . However,  for  low  void  fract- 

ions  at  flow  coefficients  greater  than  design  point  the 
single-phase  head  losses  increase  with  flow  coefficient  at 
a faster  rate  than  the  two-phase  losses,  causing  H*  to  de- 
crease with  increasing  $ . Thus,  the  randomness  of  the  cor- 

relation between  H*  and  flow  coefficient  at  low  void  fract- 
ions could  be  explainable.  For  higher  void  fractions  the  cor 
relation  shows  two  separate  trends.  For  flow  coefficients 
between  0.0  and  the  single-phase  design  point  (minimum  single 
phase  losses) , the  head-loss  ratio  increases  with  increasing 
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4>  . The  reason  for  this  is  that  the  single-phase  losses 

decrease  while  the  two-phase  losses  increase,  as  the  flow 


coefficient  increases  toward  design  point.  However,  as  the 


flow  coefficient  increases  past  design  point  the  single-phase 
losses  increase  at  a faster  rate  than  the  two-phase  losses, 
casuing  H*  to  decrease  with  increasing  flow  coefficient. 

The  flow-visualization  tests  also  led  to  some  signif- 
icant conclusions.  The  blanketing  of  the  plastic  impeller 
with  an  air  cavity  at  a 0.20  and  the  subsequent  drastic 
head  degradation  shows  that  the  flow  regime  inside  the  im- 
peller affects  pump  performance  to  a great  degree.  In 
addition,  the  flow  regime  inside  the  impeller  could  not  be 
determined  from  the  inlet  or  outlet  flow  regimes.  This  phe- 
nomenon may  well  be  the  cause  of  "unsteady"  steady-state 
data  and  inconsistent  data  readings  (reference  7-1) . It 
was  also  noted  that  for  the  same  air-  and  water-flow  rates, 
as  rpm  was  increased  (decreasing  <f>  ) the  air  cavity  took 

longer  to  form.  It  appears  that  the  high  rotational  speeds 
cause  the  inlet  bubbles  to  be  chopped  into  smaller  bubbles 
by  the  blades,  and  thereby  delay  the  formation  of  the  large 
air  cavity.  The  absence  of  the  air  cavity  in  the  bronze- 
impeller  tests  seems  to  support  this  trend  of  thought.  Un- 
like the  plastic  impeller,  the  bronze  impeller  had  blades 
which  curved  upward  into  the  impeller  inlet.  This  blade 
configuration  is  more  capable  of  chopping  up  the  inlet  flow 
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than  that  of  the  flat  plastic  impeller.  Consequently, 
throughout  the  range  of  void  fractions  tested,  the  blade 
tips  were  surrounded  by  small  air  bubbles,  decreasing  in  size 
as  the  rotational  speed  was  increased  (decreasing  flow  co- 
efficient) , rather  than  being  blanketed  by  an  air  cavrty. 

The  flow-visualization  studies  also  indicate  that  the 
location  of  the  cutwater  significantly  affects  the  flow 
pattern  in  the  volute.  During  the  two-phase  tests  using 
the  bronze  impeller  the  cutwater  created  a stagnant  water 
and  air  void,  filling  between  one-quarter  and  one-half  of 
the  exit  scroll  area.  Depending  upon  the  location  of  the 
pressure  taps  in  the  exit  scroll,  such  a flow  pattern  can 
greatly  affect  the  validity  of  the  head  readings  across  the 
pump. 

Finally,  the  importance  of  accurately  determining  the 
slip-factor  should  be  noted.  The  value  of  the  relative 
outlet  angle,  determined  from  the  slip-factor  correlation, 
is  vital  in  determining  the  theoretical  characteristic 
curves.  For  these  experiments  32  was  determined  for  a partic- 
ular operating-point  flow  coefficient  and  was  used  as  a con- 
stant throughout  the  tests.  It  is  generally  recognized  that 
32  varies  with  void  fraction  and  flow  coefficient,  although 
no  direct  correlation  is  available  at  the  time. 
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CHAPTER  8 

RECOMMENDATIONS 

A number  of  options  are  available  for  improving  upon 
the  results  of  this  experiment.  First  of  all,  I believe 
interesting  results  could  be  obtained  by  machining  off  the 
shroud  of  the  bronze  impeller  so  that  the  flow  through  the 
blade  passages  could  be  observed.  In  addition,  as  a later 
modification  the  scroll  and  cutwater  should  be  redesigned 
for  optimum  first-quadrant  operation  of  the  bronze  impeller. 
Tests  could  then  be  run  at  speeds  of  approximately  1700 
rev/min  and  for  void  fractions  of  0 to  1.0.  This  would 
provide  a more  complete  and  comprehensive  data  base  and  also 
provide  the  opportunity  to  observe  the  flow  within  the  im- 
peller passages  for  the  complete  range  of  flow  conditions. 

The  results  from  such  a program  should  be  able  to  define 
more  clearly  the  correlation  of  H*  with  void  fraction,  flow 
coefficient,  and,  possibly,  other  applicable  parameters. 
Specifically,  specific-speed  is  one  such  parameter  which 
should  be  investigated  to  determine  its  effect,  if  any,  on 
the  head-loss  ratio.  Finally,  a high-speed-strobe  photo- 
graphic study  of  the  deviation  angle,  6,  at  the  blade  tips 
could  be  done  by  injecting  dye  into  the  flow  or  attaching 
directional  threads  to  the  blade  tips.  Such  a study  could 
provide  more  information  on  the  change  in  6 with  void  fraction 
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and  flow  coefficient. 
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APPENDIX  A 

THE  BUSEMANN  SLIP  FACTOR 


The  Busemann  slip  factor,  nQ,  is  defined  as  the  ratio 
of  the  blade-tip  tangential  components  of  the  absolute 
velocity  corresponding  to  the  relative  outlet-flow  angle, 

B2,  and  the  blade  outlet  angle,  B2-  Busemann's  theory  applies 
to  two-dimensional  vanes  curved  as  logarithmic  spirals 
(blade  angle  B'  is  constant  for  all  radii).  Mathematically, 
the  Busemann  slip  factor  can  be  written  as 

oB  = (A-B<(>2tan  (90°-B£)  )/(l-*2tan  (90°-B£)>  (A-l) 

where  A and  B are  functions  of  r2/r^,  B2»  and  z.  Specific- 
ally, B = 1.0  and  is  constant  for  all  conditions  if 


2 

— > exp  ( ( 2ttcos  (90°-B  ' ) ) /Z) 
rl 

This  criterion  can  also  be  applied  to  other  than  logarith- 
mic-spiral vanes  if  B2  is  used  instead  of  B*.  Applying  this 
criterion  to  the  plastic  and  bronze  impellers  yielded  the 
following  results: 


Plastic  Impeller 


B'  = 46<> 

Z = 24 

r2  = 4 3/32  inches 


— = 1.926  > exp  [ ( 2ttcos  44°)/24] 
rl 

1.926  > 1.207 


r^  = 2 1/8  inches 


Therefore,  B = 1.0 
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e 


r 

r 


1 

2 
Z 

2 

1 


Bronze 

25° 

5 

3 21/32  inches 
1 13/32  inches 


Impeller 
r ? 

— - 2.60  >_  exp[(2ircos  65°)/5] 
rl 

2.60  > 1.70 
Therefore,  B = 1.0 


Similarly,  the  value  of  A depends  on  6^  and  z only  and  was 
found  to  equal  0.92  and  0.76  for  the  plastic  and  bronze 
impellers  respectively  (reference  A-l) . Substituting  these 
values  and  the  appropriate  flow  coefficients  into  equation 
A-l  yielded  the  following  results: 

Plastic  Impeller 


0.92  - (1. 0) (0. 5178)  tan  44° 
1 - 0.5178  tan  44° 


0.84 


Bronze  Impeller 


a = 0-76  - (1.0)  (0.03434)  tan  65°  = Q 741  = ^2 

B 1 - 0.03434  tan  65°  C' 

4>2 

From  these  ratios  and  the  information  contained  in  the  out- 
let velocity  triangles,  relative  outlet  flow  angles  of  41.76° 
and  12.5°  were  obtained  for  the  plastic  and  bronze  impellers 
respectively. 
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APPENDIX  B 

DETERMINATION  OF  THE  SCROLL  CONFIGURATION 

Since  no  performance  data  on  the  plastic  impeller  were  avail- 
able, the  first  step  in  determining  the  scroll  configuration 
was  to  assume  a (j>be*  Using  Busemann's  slip-factor  correla- 
tion, the  geometry  of  the  plastic  impeller  and  4)^  = 0.5178, 
the  relative  outlet-flow  angle,  B2,  was  found  to  be  41.76° 
(reference  1)  (see  Appendix  A) . From  the  geometry  cf  the 
blade-outlet  velocity  triangle  (see  figure  4-1) , and  employ- 
ing the  above  values  of  4>jje , $2  and  the  definition  of  flow 

C 

coefficient,  , _ _m  the  following  relations  were  arrived  at. 

4>  - u • 

C 9 

4,-  5 -iM  = 0.5178  or  C . = 0.5178u» 
v2  u_  m2  2 


r = _n»2 

z cos  39.05° 


= 0.6667u, 


Now,  the  scroll  was  designed  to  maintain  the  value  of 

C2  throughout  the  scroll.  Since  the  total-outlet  area 

through  which  Cm2  flowed  was  measured  to  be  6 square  inches, 

. 2 

the  total-outlet  area  through  which  C 2 flowed  was  (6  in  ) cos 

2 

39.05°,  or  4.66  in  . Since  the  plastic  impeller  had  24  exit 

channels,  the  area  in  each  channel  through  which  C2  flowed  was 
2 


4.66  in 
24 


= 0.1942  in  . 
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Thus,  for  C2  to  be  maintained  throughout  the  exit 

2 

scroll  the  scroll  area  must  increase  0.1942  in  every  15° 

of  arc.  Since  the  blade  and  scroll  channels  were  both 

0.5  inches  deep,  the  scroll  was  designed  to  increase 
2 

0.1942  in  _ q.388  inches  in  radius  for  every  15°  of  arc. 
0.5  in  J 
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APPENDIX  C 

CALCULATION  OF  TWO-PHASE  DYNAMIC  HEAD 
The  two-phase  dynamic  head  across  the  pump  is  equal  to 
the  sum  of  the  liquid-  and  vapor-phase  dynamic  heads. 

. rPL2CL2  J Pv2Cv2  , rPLlCLl  . PvlCvl  , 

Apdyn  tp  '-25—  + — 1 " '“Tg—  + 2g“ 


Ap. 


dEn  *P-  = AH 
^tp 


dyn  tp 


fPL2CL22  . Pv2Cv2  , 
' 2^cptP 


rPLlCLl2  . Pv^vl2, 

l ^ ~ J 


2g  p 
^c  tp 


29cptp 


where 

PtP  : 
Therefore, 


p sa  + ( 1-a) pT 
v _L 

(1-a)  + sa 


s = Jv(1-U) 


JLa 


,^L2. 2 n ,0v2 .2 

PL2(Al2}  v2  av2 

AHdyn  tp  “ f~2g~P~  * 2g  p } “ [ 


However , 


a = 


c tp 


.^Vlv 

pvl  A * 
vl 

2gcPtp 


c tp 


f°L1i 

L1 V 

2g  p. 

^c  tp 
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Resulting  in 

AH  . 

dyn  tp 


„ , L2  ,_v2_. 

PL2  (l-«)  PV2(°AT2) 

2gc°tp  2gcptp 


PL1( (1-«)At1)  ( Pvl(aATl) 

2g  7 P + 2g  p . ^ 

tp  tp 


In  performing  the  two-phase  dynamic-head  calculations 


it  was  discovered  that  the  contribution  of  the  vapor-phase 
was  negligible,  and  was,  therefore,  neglected  for  subsequent 
calculations . 


I 
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APPENDIX  H 


THE  DRIFT-FLUX  MODEL 


The  drift-flux  model,  developed  by  Zuber  and  Wallis, 
satisfactorily  accounts  for  the  influcence  of  mass  velocity 
on  the  void  fraction  and  is  useful  in  the  bubbly-,  slug-, 
and  churn-flow  regimes  (reference  3-2) . The  details  of  the 
drift-flux  model  are  given  below. 

J 


a = (- 


1 . 2 ( JT  +J  ) + 0.35 (gd) 

L V 


172 


+ = upflow 


- = downflow 


(D-l) 


_ *v 


J = superficial  vapor  velocity  = . - = 


V T 


*L  °L 

IT  = superficial  liquid  velocity  = - . ■ - = j— 


L T 


d = inlet-pipe  diameter 

Substituting  for  J and  JT , equation  D-l  becomes 

V ii 


a = [— i ITT- 1 

1 . 2 (Q  +Qt ) - (0.35) (gd) 


(D-2) 


*V  “L'  y - • -*T 

After  some  mathematical  manipulation  equation  C-2  can  be 
transformed  to 


Qv  “ 


1 . 2uQT  - 0 . 35 (gd) 1/2A_ 

(1  - 1.2a) 


(D-3) 


A 


100 

Thus,  for  a desired  void  fraction  and  water-flow  rate  an 
appropriate  air-flow  rate  was  calculated  using  equation  D-3. 


t 
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APPENDIX  E 


COMPUTER  program: 


CALCULATION  OF  TWO-PHAGE  DYNAMIC  HEAD 

0004.1 

D E N K A T = 6?. 4 

C00CR 

A R e,  A = 0.3541 

00 1 4 R 

HO  = 51 

001CS 

Cr  = 531.44 

0024R 

HRITE  (5,4) 

00  38  a 

8 

FOHNAT  (141,  5 X, * TTST  ho.  •,  5X,  • 

0038P 

1 *,5X,  * VOID  F .< ACTIO"'  *,  SX,  K>M 

O038H 

2DENSITY  * ) 

0DBER 

DO  5 I = 1,HD 

O0C6R 

BEAD  (8,6)  HUH,  QKAI,  OAIP,  S,  tp'T';,  VOID 

OlOAfl 

6 

FORNAT  (110,  2F12.5,  T10,  2P12.4) 

0128R 

QTnr  = OH  AT  ♦ (.AIR 

0 1 34  h 

FC05FF  = (l;TOT/FLOAT(N))  * CF 

014C8 

DEN  AIR  = (PRESS+14.7)  *144.  /(53.3«  » 53C.) 

01708 

SPA I R = C’AIP  / AREA 

017CS 

SPWAT  = 0 H A T / APF» 

01888 

SLIP  = STAIR  » (1.-V0IP)  / (SPWAT  * V'lT) 

OUCH 

TPDEN  = ( DEN  AIR* VOID *SII P ♦ (l.-VOT'i)  * 0EN«*T) 

0 1 ACS 

1 /(I. -VOID  ♦ SLIP  * VOID) 

01F4R 

DYNHD  = (1076. 5 305*  Oil  AT*  • 2 )/(  TPDE  X *(  1. -VOID)** 

02288 

HRTTF.  (6,10)  HUP,  FCOE'F^,  VOTD,  Of HHf,  T FT"* 

0264R 

10 

FOHNAT  ( 5X,  I ID , 4(F17.5,5X),  /) 

0280R 

5 

CONTINUE 

07101 

END 

'X' 


CDS FFrCT- 
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SINGLE- PHASE  DATA  CURVE-EIT 

030«h  omt'isi'"'  hot 

0304ft  DIIENSION  FLCu£F()'>),  HSP<35>,  CH  F HO),  PEI*  Z'.H  ( )0 ) 

0304S  DH“)iEIOft  C(TFF(7) 

0004ft  <13  = 2) 

030CR  AND  = 2'». 

03 1 4 K NCtF  * 7 

03 ICR  DO  5 1=  1,  ID 

0024ft  REAP  <8,6)  FlC<"«F(I),»'Sf  (T) 

0078  3 PLOT  (1,1)  * )LCO‘:r(l) 

03AAR  PLOT  (2,1)  = !S( (T) 

03DCb  6 FOHMT  (2F12.b) 

03EA3  5 COST  INUF 

OOFAft  CALL  L3FIT  <4!),  HOFF,  FLCOFF,  HSR , COEFF) 

oio ak  wniT2  (5,))  co’.rr 

0128ft  9 FOR*  IT  (//,  • THP  COF.FFICI  -NTS  ARE  : •,  7K17.5) 

0158R  WRITS  (5,15) 

016CR  15  FOftflAT  (////,  5X,  * * , 2U,  *HSF*  , 17X,  ‘CHECK  * , 21X,  ‘PERCENT*  J 

01A6R  C 

01A6S  TCTE’Rh  * 0. 

01AER  C 

01 AEB  DO  10  J = 1,  ND 

0 1 B6h  CHECK(J)  * COEFFd)  ♦ C0EFF(2)*FLCCFF(J)  ♦ COEfFO) 

01B6R  1 •FLCOEF(J>**2  ♦ COEFF(4>  • FLC0EF(J)»«3  ♦ C0EFF(5)*FLC7EF(  J) 

01B6R  2 **4  ♦ COEFF(6)»FLC3EV( J)*»S  ♦ COEFFO)  * ELCO~F(J)**6 

02CAP  PLOT  ( 3, J ) = CHEC<( J ) 

02FCR  PERCEH ( J ) « CHECK(J)  / HSP(J) 

0338R  ERR  * 1.-  ABS  (PFRCEN(J)) 

0364ft  TOTESB  = TOTEPF.  ♦ ERR 

0370ft  WRITE  (5,20)  FLCOEF(J),  HSMJ),  CHECK(J),  PrRC;  N ( J ) , ERF 

040CK  20  FOR1AT  (/,  E13.5,  13X,  E13.5,  10X,  E13.5,  2DX,M9.u,  15<,F1".4) 

0444R  10  CONTINUE 

0454R  C 

0454  R AVERK  = TOTERR/AND 

0460R  WRITE  (5,25)  TOTERF,  AVEKR 

0484R  25  FORPAT  (////,  * THE  TOTAL  ERROR  IS  *,F10.4,‘  ANi  THE  A V 7 . 

0484R  1 ERROR  IS  *,  F10.4) 

04DAR  CALL  QPICTR  (PLOT,  3,  29,  31(2,3),  3X(1),  3LAPEK -1 13 3 ) ) 

051ER  PAUSE 

0524R  CALL  jPICTP  (PLOT, 3, 29,  01(2),  QX ( 1 ) , 9LABFI,  ( -1103)) 

3566  R PAUSE 

056CR  CALL  OHCTK  (PLOT, 3,  29,  QK(3),  QX(1),  3LAREU -1 10 3 ) ) 

05AER  PAUSE 

05B4P  END 
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SINGLE 

00043 
003  43 

00  34?. 
0034r. 
0034? 
0034? 
0034  R 
0034? 
00343 
0034? 
0034? 
J034F 
0034H 
03343 
03603 
009CR 
00*43 
03*43 
03*43 
03*C? 
30E2R 
03E28 
01063 
010ES 
31 36h 
&136P 
0136.3 
313EH 
01463 
01683 
018*3 
01928 

01  £83 
01F43 
32203 
02308 


PHASE  DATA  CURVE-FIT  (continued) 

SUHnOUTJSK  I.SF  I T(  N .TS  , NPAKA1,  ,•  , Y , (»»»*«) 

C 

y,  r»?\p,  i,  xpow?.*,  xr,  n,  yk 

1 9 TEG  HP  N'PTG  , PAPA",  I.  K,  l**?**,  IrKP 
C 

DIMENSION  X(NPTS),  Y(‘PTS),  PARA*(NPARM) 

C 

C TO  USE  A VN.VJE  Or  SPAR  AH  IprAT-  P THAN  10,  TH~  IIC-s-  f H " U LO  CHAN  ;E 
C THE  LENGTH  Or  * TO  AT  Lf  *3T  VPARA*»*2#  AND  THE  LF”GT«  OF  X °0i<  k 
C TO  AT  LEAST  2*<NP r RAK-1 ) . 

DIMENSION  A(IOO),  XPOUEM  16) 

C 

C CHECK  FOP  ARGUMENT  ERRORS. 

If  <fiPTS  .SE.  ’.'PARA1'  .AND.  M P A P A N .CT.  0)  GO  TO  ID 
WRITE  (5,  1001)  NPT3,  NPRRAY 
RETURN 
C 

C ZERO  ARRAYS  PEFCRr  SiJv*IN5. 

10  1)0  20  I = 1,  NPAPAY 
20  PAM»d)  = O.OE'I 
C 

IXPNAX  = NAKOO,  2 * (NPARAX  - 1)) 

DO  30  I = 1,  TXP"AX 
30  TrOWffcU)  = O.OFO 

c 

C CONFUTE  SUNS  OF  POWERS  OF  X «Nf  OF  POWERS  OF  X TT"eS  T. 

DO  SO  K = 1,  NETS 

XP  = 1.0E0 
XK  = X ( K ) 

YK  - Y ( K ) 

DO  40  I = 1.  IXPrAX 

IF  (T  .Lr'.  NP'.HAT)  PKiAKd)  = PASAPd)  ♦ X P * Y <C 
X P = XP  * <K 

XPOWES(I)  = XPOWE3  ( I ) «■  XP 
40  CONTI  NiJi 

so  oxrr  u*' 
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SINGLE- PHASE  DATA  CURVE-PIT  (conlinued) 


D244rt 

C 

024  4 m 

C 

coKi'iirt  Ot-tKirir.! k r.j-'U  ■< r,->  co-  mri 

0244M 

C 

HAS 

; a a i>  str'IiTOR’  . 

0244m 

a.  ( i > = arr, 

026  AM 

if  (vp.r.s'  ,r o.  i)  so  ro  9r 

0270k 

C 

0270K 

DJ  RO  I = 2,  Vf'ARA« 

0279K 

X.»  = XKOMEI  ( T - 1) 

0296  R 

"J  60  * = 1,  I 

029EK 

C 

029E.it 

C 

A(K,  I + 1 - A)  = XP 

C29ER 

JK  = ( K+’N'RAHA*  * ( I -K ) ) 

02B6K 

A ( JK  ) = XE 

G2CER 

60 

CONTINUC 

02DEK 

C 

02DEh 

JL  =<NRAftAN  ♦ I -2) 

02F2B 

XR  = KP3WER  <JI.) 

030AR 

DO  70  X = X , NFARAf! 

0312R 

C 

0312R 

C 

A( NPARAM  ♦ I - X,  K)  * XP 

0312B 

JH  =(I  ♦ (NPAM»  -1  ) * n 

0 32  AR 

A ( JH ) = X° 

0342R 

70 

CONTINUE 

0 356  ft 

80 

CONTINUE 

036AS 

C 

0 36  Aft 

c 

SOLVE  NORMAL  E3UATT0SS. 

036  AR 

90 

CALL  SINC(A,  r ARAN , NPAPA*,  IEFR) 

0390R 

IK  ( IEHR  • NE • 0)  WRITE  (5,  1002) 

03B2R 

RETURN 

03BAM 

c 

03BAR 

1001 

FORKXT(?3u  LSFTT : ARGUMENT  ERROR,  •’111) 

03E2R 

1002 

FORMAT!  392  LSFITs  N'OhRAL  EOUATIONT  ARE  SINGULAR.) 

0414R 

ESC' 

- 
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CALCULATION  OF  ^ ' 


0304R 
0004H 
0004H 
0004K 
0004H 
C314K 
001CR 
0324  k 
00303 
0338r 
0336  k 
004CR 
004CE 
009CR 
009CH 
OOBOK 
0QB03 
OOBOS 
013CH 
313CR 
0150R 
OlbOR 
01C8R 
01C8R 
01D0R 
02A4R 


rtPth  ' 


FOR  THE  PLASTIC  IMPELLER  DATA 


Nj-.doci,  •...  r o t < i w,  o » , \ ( 1 3 o ) > r-s  snoro,  vu  m 3) 
DI/EN3IDA  'JLID(100),  !::.*<  * T ! ( 1 00 ) , MV'!)),  n"I  ( - :,r.)  , '»  r ■ ,n-<  vf.) 
D I NX  IS  5 1 ON  ,Mi(1i)C),  ;.rfJT(1f>0) 

AREA  - C.2j»24.  /1U4. 

AA  = 1.0 
N3=)1 

R=-1.1?0ri14 

CF-671.75 


WRITE  < b , b ) 

F 3 i;  1 A T ( 1 k 1 , / / , 3 1 < , • C ?.  L C l*  L A i N T 
1RUBSL?  PU'P.') 


TH E TWO-PHASE  F!7V  HINTTTO" 


10 


writei*. 9) 

FOR  N AT  ( ///  » 4 X * * T~ST  ’10.  1 , * TOTAL  PLOW,  RV,  • *;  r.' }•",  • , 1 ;;  , 
1 * V 3 1 L>  • , 12*,,S^P,<3X,,'t'P^•^FICIAL,,  b)i  , • SUPF'. 1 I n 
2 • TWO-PHASE  * , 3X , • TSF0PET1CA1 ■ ) 

WRITE  (b  ,13) 

FORN  AT(  IbX  , *CU  .FT/SF-C  * » 7X  ,•  PP*  • ,3X,  • FRACTION  • ,1 1 k A » T*V  , 

1 *»iL  AIK’.BX.’VEL  LID.*  ,10X  , * FUNCTION  • ,27  , • TWO-PK  AM  S ” :.D  * ) 

no  b i = i,  no 

READ  (8,4)  NUN(I),  0»AT(T),  0*  IS  ( I ) , N ( I ) , FEES'  (I),  V3IJIT) 
FOR*! AT  (TIP,  2F12 • b , Til,  ?H2.4) 


02C2rt 

02C2R 
02C2R 
02CAR 
0 30  ER 
030ER 
033A3 
0356R 
0362K 
036ER 
03AAR 
03DAR 
03FER 
04  52  R 
04AER 
04AER 
05S6R 
ObCAR 
ObCAh 
ObCAK 
30CEH 
06FAR 
070AR 


C 

C CALCULATE  THE  DENSITIES  OF  WATER  ASP  AIR. 

DEN  WAT  = 62.4 

DENI.  I r!  (I  ) = (l-p.KSS(I)  * 14.7)  *144.  / (5  3 . *b  '3  . ) 

C 

AIRAA3*  'JAlR(l)  * JLNaIP(I) 

U ATI  AS  = 3WAT(T)  * DEN » AT 
TOTN  AS  = A I fi  *A  S ♦ K ' r*AS 
X = AIHNAS  / tot-;.:: 

QTOT ( T ) = 3 A I k ( 7 ) +3WAT(I) 

spair  = Tori  a?;  * x / (DEVAir  (i>  ♦ ape.a) 

SPWAT  = T3r«AS  * < 1 . -r ) / (PEN'WAT  * AREA) 

SUPd)  = SPA1F  * (l.-VOTP(I))  /( SPK  AT  * VOTD(T)) 

A(  I ) = (VOIPd)  /(1.-V0IIHL)))  * DM'M’U)/  DF'W'T 

PHI  C I ) = ( 1.  + KI  ) ) * (1.*  “T.TP(I)**?  » A ( T ) ) / ( 1.  * SLTP(T)*A(T)) 

1 •»? 

KTMTP(T)  = AA  ♦ » PHMT)  * (0TOT(T)  / (FLOAT(*(T))))*SP 

C 

WHIT1'  (',,7)  JHM  I ) .jTOr(I),N(I),V'>ri)(T),'U  U-(T)  ,~.l  M3, 

1 SPWAT,  PHI(I),  HT'irPU) 

7 FUH'AT  (/,4  y,  Vl,  (.\,  FlO.b,  110, 4X,  f,  { FI  2 .5,  3 X ) ) 

5 70NT I NIJ1’ 

END 
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CALCULATION  OF  II*  FOR  THE  BRONZE  IMPELLER  DATA 

30  34  , LT1  (7,  :>1 ) 


0004* 

OK-  v .fO1'  r LOT  • 

(-,=  1 ) 

0304  ( 

■)i  i.ivstov  :<oc: 

(4) 

0004.-. 

DATA  SCI  /O.O, 

0 f ) • / 

001  bH 

c 

00164 

;n  ii  ti  = i,s 

0 0 2 0 R 

00  12  JJ  = 1 ,S1 

0023K 

PLOT?  (II, IJ)  = 

?.7? 

004AR 

12  CONTINUE 

035AP 

11  CONTINUE 

036  A H 

A A = 1.0 

0372R 

AREA  = 0.0681 

007  AS 

RD*!3  1 

0082. < 

B=  -4.611 

008E.R 

CF 1 = 639.44 

0096  R 

CF2  = 26.36.774 

009ER 

ISCL  = -2 

O0AAR 
00AAP 
ODbEh 
03BER 
03RES 
01 32R 
01  32R 
013AR 
0 1 8 6 R 
01A4R 
01ACH 

0 1D0R 
01DCR 

01  OCR 
01E8R 
01F4L 
0218R 
0230. < 
0278  R 
0290R 
0290R 
02BCK 
02DCR 
02DCR 
0370R 
0390R 
0 3 BOR 
03B0R 
03D2R 
0 3 F 4 f . 
04 16  K 
04  38  P 
046AK 
047CR 
049EH 
04CEP 


WRITE  (S,?3) 

23  F0«?1  AT  ( 1H  1 , 37,Ttor  S3.*,  2X,  'VOID  FRACTT0V,  M, 

1 • FLOW  COEFF*.  6.<  , 'HTHTP  ' , 12X,  *HTP',  10\,  • H^PSP  1 , 

2 12X,  'HSP*,  13X,'R»TTn' ) 

C 

DO  & I = 1,  40 

READ  (P,£)  S Uf!  , ?AI9,  9,  PRESS,  VOID  , Dt.LP 

6 FORMAT  (113,  2F12.r>,  110,  3F12.S) 

D E 4 W A T = 62.4 

DK4AIR  = (PRESS  + 1m .7  ) *144.  /(S3.34  * S3  0.) 

OTOT=  gAIR  ♦ CWAT 
C 

SPAIP  = 2 A I P / AREA 
SPWAT  * 2WAT  / ARP? 

SI.IP  = SPAIR  *<1.  - VOIP)  / (SPWAT*  VOID) 

A=  (VOID/  1.  - VOID)  • DF4A.IB  / DEKWAT 
PHI  = (1  .♦?)•(  1.+SLIP**2*A)/(  1.+SUP*A  )*»2 
FLCEF=  (2TOT/FLOAT(!i))  *CF1 
C 

HTHTP  --  AA  + H»PHI  «< QTOT/FLDATI 4 ) ) *CF1 
HTHSP  = A A ♦ B*OT0I/FL0AT(N))  »CF1 

rtSP  = .41315  ♦ .673  6 1 *FLCEF  - 5 .77  19»FLCEF«*2-  6 14 ,27*FLCEF**3  + 
17931. 1*FLCEF*«4  -46S2f>.  *FLC£F*»S  ♦ 106 8 1 3 . « FLCE)  * »f, 

HTP=  (DELH  /FL0AT(4)  **?)  *rF2 
RATIO  = ( HT.ITP-HTP )/( HTHSP-FSP  ) 

C 

PLOT  (1,1)  = FLSEF 
PLOT  0,1)  = HT'ITP 
PLOT( 3, I ) = ir1 
P LOT (4,1)  * 'ITTSP 
PLOT ( S ,1 )=  MSP 
PL0T(6,I ) * RATIO 
PLOT ( 7, I ) * VOIO 

IF  ( FLCEF  . Lk.C.is)  i>r.OT?(  1,T)«RATTP 

1F(FICEF  ,2T.C.C'8.AND.  FLCPF  .IF. P.11)  PLOT:  ( 7 , 1 ) = “ » TIC 


T 
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CALCULATION  OF  II*  FOR  THE  BRONZE  IMPELLER  DATA  (continued) 

0S0CR  IFtFLCtF.ST.O.  11  .MI'.FLCFF.LE.O.IU)  FLOT2  ( ? , I ) = F 'TTC 

0S4AK  IF(FLCEF.3T.0.14)  fLOT2(U,I)=  PATIO 

0S7AH  PLOT2(S,I  ) =V 01 C 

0S9CR  C 

059CR  WRITE  (E,18)  NU 1 > VOID,  FLCKF ,HTHTP , HTP.  HTHSf  , IT?,  'ATM 

OSFOP  18  FOHKAT  (3X,  110,  7 ( F10 . S , SX ) , / ) 

OSOCR  5 CONTINUE 

061CR  CALL  CPICTR  (PLOT,  7,  «D,  CT<  2,3,4,'.),  3X(1),  .>LAPI  L< -1 70  7 ) ) 

0664R  PAUSE 

066 AR  C 

066 AR  CALL  CPICTR  < PL0T2 ,S , NO.  OY ( 1 ) , 3X(5),  3LAPFL( - 100  7 ) , 

066 AR  10ISCU-2)  ,0XSCL(XSCL)  ) 

06  DU  R PAUSE 

06DAR  CALL  OPICTR  (PL3T2,5,ND,  0H2),  OX  f S ) , 3LAREL( -1 00 ?) , 

06DAR  1QISCL(-2),0XSCL(XCCL) ) 

07UUR  PAUSE 

074AR  CALL  CPICTR  (PI0T2,5,RP,  OKI),  QX(S),  3LARKL< -1 00.0  , 

074 AR  1QISCL(-2),3XSCL<XSCL) > 

07B4R  PAUSE 

07BAR  CALL  QPICTP  ( PI 0T2,5 , ND,  CT(4),  3X(S).  JLARFL( -1  DO?) , 

076AR  1QISCL(-2),0XSCL(XSCL)> 

0B24R  PAUSE 

08  2 AR  END 


t 
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